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SYNOPSIS 


During recent times effort has been made to understand the pairing mechanism in high-T c 
superconductors. Some theorists believe that an effective electron-electron attractive interac- 
tion mediated by spin fluctuation is the possible source of pairing mechanism[l], Two-magnon 
Hainan scattering is a very useful tool to probe the low-energy excitations in magnetic solids 
and hence can provide information regarding the nature of pairing correlations. 

In the sixties extensive studies have been carried out to understand inelastic light scat- 
tering phenomenon in antiferromagnets such as RbMnF 3 , MnF 2 , K 2 NiF 4 etc. [2,3] and a good 
theoretical understanding exists for these systems with spin S > 1 , within the interacting 
magnon theory for the Quantum Heisenberg antiferromagnets (QHAF). Even though the 
two-magnon Raman scattering has been studied extensively within the QHAF, there have 
been so far only limited attempts to systematically study this within the Hubbard model. In 
spite of the recent success by Shastry and Shraiman[4] in obtaining the Fleury-Loudon (FL) 
Hamiltonian by performing a strong coupling expansion in powers of t/U, t and U being the 
hopping term and the on-site Coulomb potential, so far no systematic expansion (a formally 
weak coupling one) has been carried out which can continuously interpolate between the 
small-/* and larged* limits. 

Recently the study of two-magnon Raman scattering gained interest since the experi- 
ments performed on La-H'uO* [a] (a spin -5 Heisenberg antiferromagnetic system in two di- 
mensions) reveal a two-magnon Raman scattering linewidth much too broad to be explained 
within the classical theory. Moreover scattering intensities are observed in geometries, which 
in tlu* classical approximation, should yield no scattering. The third issue, which seemed 
puzzling was that tin* two-magnon spectral weight persists all the way upto 6.7 whereas one 
would expect the spectrum should have a cut-off at 4.7 because this is the maximum possible 
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energy that the non-interacting magnons can have. 

Some of the discrepancies between the theory and the experiments mentioned above have 
been clarified to a certain extent by Singh et al.[ 5] by incorporating quantum spin fluctuation 
which is inherent to low-spin (S = 1/2) and low-dimensional (2D) systems. Recently it was 
emphasized by Nori et a/.[6] that it is necessary to invoke additional processes where they 
took into account magnon-phonon interaction which can contribute to both the two-magnon 
Raman linewidth and the finite intensity in the classically forbidden channel. They supported 
their arguments by computing the Raman cross-section using a S — 1/2 Heisenberg model 
with random static exchange integrals on finite size lattices using Exact Diagonalization (ED) 
and Quantum Monte Carlo (QMC) techniques. We have pointed out that these calculations 
do not sufficiently explain the explicit effects of the disorder upon the high energy magnon 
modes with energy ~ 2 J which are mainly responsible for the Raman scattering. 

Very recently Chubukov et a/.[7] have studied extensively the two-magnon scattering 
process in the resonant regime where the incident photon frequency is of the same order as 
the band gap in cuprates (both are ~ 3 ev). They found a strong divergence of the Raman 
vertex close to resonance which was useful in explaining the experiments done on S^CuChCh 
and YBa-V 'uriOr.-rx [8], 

Apart from the exchange disorder which presumably is arising out of zero-point lattice 
fluctuations, there are many other correlated electron systems which are intrinsically dis- 
ordered. The interplay of disorder and spin-fluctuation effects is disordered antiferromagnets 
have been of intense interest in recent times[9]. In view of that we carry out a comprehens- 
ive study on two-magnon Raman scattering in disordered antiferromagnets which yields a 
detailed information about the spin fluctuation spectrum in presence of on-site disorder. 

In chapter I we introduce the topic and provide the necessary background. It contains 
an overview of the theoretical and experimental status of the field. 
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In Chapter 2 we develop a systematic diagrammatic scheme to study two-magnon 
Raman scattering in the Hubbard model. There the fermionic structure of the magnon in- 
teraction vertex is studied at order l/Af where A f denotes the orbital degeneracy. The full 
two-magnon propagator is then evaluated by including the interaction between the magnons 
at the ladder sum level. The scattering intensity is then obtained from the imaginary part of 
the propagator. Finally the calculations are extended to the intermediate coupling regime in 
order to make connections with the experiments done on cuprate insulators. An attempt has 
been made to incorporate zero temperature magnon damping which yields a better agreement 
with the experiments performed on La 2 Cu 04 - 

In Chapter 3 we study the effect of hopping disorder in a Hubbard Antiferromagnet us- 
ing both perturbation theory and exact eigenstate method. Here the lattice distortions arising 
out of quantum and thermal fluctuations are considered within adiabatic approximation (or 
Born-Oppenheimer approximation) which neglects the fluctuation of the phonon field. It is 
shown that while the low-energy long wavelength modes are only weakly affected, the high 
energy modes are strongly renormalised because of a cooperative effect arising from local 
correlations in hopping disorder. We also discuss the implication of this highly asymmetric 
magnon energy renormalization for two-magnon Raman scattering in cuprate antiferromag- 
nets. 


Chapter 4 deals with obtaining the Raman scattering intensity numerically in presence 
of hopping disorder. A consistent scheme of mode assignements is used which essentially 
establishes a one-to-one correspondence between the pure and disordered system and we be- 
lieve it to be valid in the weak disorder limit. It is shown that even for a substantially weak 
disorder, several features like the asymmetry with respect to the two-magnon peak, linewidth 
to peak position ratio and the spectral weight persisting all the way upto 6 J agree well with 
the experimental data. 


in Chapter 5 we have calculated the Raman scattering intensity in an antiferromagnet 
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in presence of random on-site energies. It was found[9] that the magnon modes ace damp-ed 
due to scattering of magnons by disorder. This damping was calculated by considering the 
vertex corrections and shown to be proportional to the fourth power of the disorder strength. 
We obtain the magnon damping in a self-consistent manner and use this to calculate the 
Raman lineshape for various disorder strengths. 

In Chapter 6 we develop a real-space formalism to study the two-magnon propagator 
which allows the calculation of the Raman intensity for disordered systems involving hopping 
disorder, potential disorder, impurities, vacancies etc. 
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Chapter 1 
Introduction 


1.1 General Discussion 

Experimental observation of the inelastic light scattering in liquids and solids has a long 
history dated since 19'2S due to Raman [1, 2] which won him the Nobel prize in Physics. 
He observed that when the incident light of frequency u interacts with the system, then 
the scattered light has a frequency w 3 , which is either greater (Stoke’s Scattering) or lesser 
(anti-Sfoke\s scattering) than u Also there can be a possibility where the incident light 
suffers no frequency shift due to scattering and is called Rayleigh scattering. Brillouin[3] 
and Man<lelstam[4] predicted similar scattering effects from sound waves in dense media. 
The term ’Raman scattering’ is now used to describe light scattering from various kinds of 
excitation like phonons, magnons etc. In my thesis, I have concentrated only upon light 
scattering from magnons or spin waves which are the lowest-lying collective excitations that 
occur iii magnetically ordered material. 

1.2 Two-magnon Raman Scattering from Pure Antifer- 
romagnets 

It was shown by Elliot and loudon[5] that the mechanism responsible for light scattering from 
magnons is the electric dipole interaction. The magnons in the magnetic material create a 
spatially periodic modulation of the permittivity and the light gets scattered by the fluctuation 
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in permittivity. The first observation of Raman scattering from magnons was carried out by 
Fleury[6] where he studied the temperature-dependent light-scattering spectra from one and 
two-magnon excitations in antiferromagnetic FeF 2 - Since the incident photon has essentially 
zero wavevector, one-magnon scattering arises from magnons with wavevector near the cen- 
ter of the Brillouin zone (BZ), while for two-magnon excitation, the magnons have equal 
and opposite momenta. Soon many experiments followed to observe Raman scattering from 
magnons in other antiferromagnets like MnF 2 [7], RbMnF 3 [8] and NiF 2 [9]. An apparently 
surprising result emerged from these experiments which is the two-magnon cross-section is 
usually larger than the one-magnon cross section. This can be trivially explained by noting 
that the scattering cross-section is proportional to the magnon density of states (DOS) and 
the DOS is very small for magnons with wavevector close to the center of the BZ. I will be 
mostly interested in looking at the two-magnon scattering phenomenon in antiferromagnets. 

Elliot and Thorpe[l I] have shown that in an antiferromagnet the spectrum of the scattered 
light as a function of energy transfer from the photon to two magnons is very sensitive to 
the interaction between the two magnons. The reason is that the two magnons are created 
at two neighbouring sites and hence the interaction is rather important (Fig. (1.1)). The 
necessity of invoking the magnon-magnon interaction to explain the complete lineshape was 
experimentally confirmed by Fleury[8] on RbMnF 3 , Thorpe[12] on MnF 2 , Parkinson[13] and 
Fleury et a/. [14] on two dimensional antiferromagnets such as K 2 NiF 4 . The essential feature 
of the scattering spectrum borne out of experiments is that the two-magnon peak occurs at a 
transferred photon frequency of J(2SZ — 1), whereas the joint magnon DOS peaks at 2 JSZ. 
It indicates that the presence of the first spin deviation leads to a reduction in energy for the 
second spin deviation. 


At zero temperature, the two-magnon Raman cross-section was derived by Elliot and 
Thorpe[ll] using a Green’s function theory. They employed a refined decoupling procedure 
to linearize the Green’s function equations in a controlled manner. It is appropriate here to 
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Figure 1.1: Two-m.u'.iit mi Hainan lineshape for Rl^MnF.i without maguon-magnon interaction 
(curve A) and with mag non-magnon interaction (curve B). (From Ref. [13]) 

briefly describe the t inn try due to Elliot and Thorpe(15]. 

To calculate the two inngnon Raman spectrum one needs to calculate the imaginary part 
of the two-magnon propagator given by, 

G(S, 5') = (Pj Pj,) 

Whei'e Pj is the tram. verm* spin-pair operator given by, 

P I = 2 12( S r S r+ ^r Sj + s) 

and S, S' are unit verttus connecting nearest neighbours. When one writes down the equa- 
tion of motion for U\ r l f r ) following Zubarev[16], it turns out that one has to evaluate 
the commutator \lj, H H i where II U is the light scattering Hamiltonian given by Floury and 
Loudon[17], 
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Hr = A^(E inc i)(E sc i)S r .S r+l (1.1) 

r 

where Ei nc and E sc are the incident and scattered Electric field vectors of the photon. The 
commutator leads to terms like where either all of the site indices ri, r2 and r 3 

either correspond to different sites or at most two sites are same. The above product of 
spin operators must be decoupled in order to obtain a set of linear equations that may be 
solved to yield ( Pj Pj , ) . A simple minded decoupling is to replace by ±5 depending 
upon the sublattice. But this yields the non-interacting spin-wave result which has got very 
poor agreement with experiments. The deficit was taken care of by introducing an improved 
decoupling procedure. If r l5 r 2 and r 3 are all different then replace S* x by ±S appropriate 
to the sublattice, while if 17 = r 2 then replace 

s;x> ss+ 

and, 

for spins in the t-spin sublattice. Similar replacement for spins on the .(.-spin sublattice (by 
changing S to — 5 ) and also for the case ri = r 3 . This decoupling straightforwardly yields 

the two-magnon cross-section to be proportional to, 

T Gp 

Im ll + JGo. 

where J is the self-energy arising from the magnon-magnon interaction and Go is the non- 
interacting two-magnon propagator given by, 

Go(^) = yv ^ (ft 2 -4ft*) 

where <p q is the symmetry factor appropriate to the lattice and Cl g is the spin-wave energy 
obtained at the RPA level given by Q q = ZJS^j 1 - 7, 2 and 7, = \ £* e i(,:& where Z is 
the number of neighbours. The spectrum thus obtained has got excellent agreement with 
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experiments. A generalization of the calculation to finite temperature was done by Solyom[4] 
by taking into account the temperature dependence coming from the Bose factors describing 
the occupation number of magnons. He finds that the two-magnon peak shifts with increasing 
temperatures towards low energies and the peak broadens. Also it is observed that the peak 
does not disappear at the Neel temperature T}v and even continues to persist upto as high as 
47 ; v- This can be qualitatively understood by arguing that the short range spin correlation 
plays an important role in the two-magnon Raman scattering and that does not vanish at T,v. 
There were other theoretical efforts to extend the theory at finite temperature during the early 
seventies. These include a diagrammatic scheme using the Dyson-Maleev transformation[19] 
for T < -Tv [20] and the equation-of-motion approach using a modified decoupling procedure 
for T > Tv [21]. 

1.3 Two-magnon Raman Scattering from Impurity Doped 
Antiferromagnets 

The introduction of low concentration of substitutional impurity into an antiferromagnetic 
insulator modifies the spectrum of magnetic excitations. In general, the spectrum consists 
of a set of impurity modes associated with the defect and it’s immediate neighbours. The 
impurity modes occur outside the band of host excitations and are localized in the neigh- 
bourhood of the defect. 


1.3.1 Non-magnetic impurity 

Raman scattering from pairs of magnons in antiferromagnets in presence of non-magnetic 
impurities like Mg, Zn etc. were studied experimentally for both low and high impurity 
concentrations. It was observed that the position of the peak shifts towards lower frequency 
and the spectrum broadens with increase in impurity concentration. Fleury et al. [22] studied 
Raman scattering in doped antiferromagnets like KNij_ x Mg x F 3 upto a concentration x as 
high as 0.75. Also experiments were done on Zn doped systems such as Mni_ x Zn x F2[23, 
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24] and it displayed a similar concentration dependence as observed in KNi 1 _ x Mg x F 3 . The 
two-magnon peak was observable upto x ~ 0.7 but was very weak. Again the peak fre- 
quency decreased linearly with x, but the linearity persisted only upto x ~ 0.3. Some other 
groups[25] have looked at Raman scattering from mixed crystals like Fei_ x Zn x F 2 where the 
peak frequency showed a linear fall off with x upto x — 0.5 and then it roughly remained 
constant upto about x = 0.75. 


1.3.2 Magnetic Impurity 

Broadly we can classify the experimental work done on Raman scattering from magnon pairs 
in antiferromagnets diluted by magnetic impurity into two types. The first is the scattering 
from point defects i.e. at very low impurity concentration (e.g. ~ 1%.). Secondly, at high 
concentration the scattering involves excitation of magnons in each sublattice A and B giving 
rise to three kinds of pair excitations, viz. A — A, A — B and B — B. 

Scattering from point defect type had been observed in several mixed antiferromagnets 
viz. Ni doped RbMnF3[26, 27], and KMnF 3 [28] and Co in RbMnF 3 [29]. The general obser- 
vation is that two-magnon peak moves towards higher frequency with increasing impurity 
concentration. With higher concentration the two-magnon Raman spectrum exhibits three 
distinct bands. The most detailed experimental account was provided by Lockwood et a/. [30] 
where they studied the light scattering from KNi x Mni_ x F 3 for four different impurity concen- 
tration such as, x = 0.84, 0.73, 0.54 and 0.15. The spectrum showed clearly two-magnon 
bands corresponding to the creation of magnons on neighbouring pairs of Ni ions, neigh- 
bouring Ni and Mn ions and neighbouring Mn ions. The frequencies at which these bands 
occurred had got good agreement with theoretical calculations done using Ising Cluster model 
along with an improvement done which incorporated the modifications in the exchange con- 
stants with lattice parameter due to doping. They had also performed a systematic study to 
look at the temperature dependence of these bands. The bands became broader with increase 
in temperature and shifted towards lower frequencies. Qualitatively the behaviour remained 
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the same for all the bands. 

Before I close the discussion on light scattering from impure magnets, a very brief discus- 
sion of the Ising Cluster model is presented here. It provides a good description of frequency 
dependence of the scattering process in pure materials essentially because the zone boundary 
frequency of the magnons in antiferromagnets is given quite precisely by the Ising model. 
These zone boundary modes dominate the two-magnon scattering process because of their 
large DOS. Consequently we expect it to provide a substantial information of the two-magnon 
Raman scattering in these systems. 

For high impurity concentration, the Raman intensities of the various bands in the spec- 
tra reflect the probability of finding the corresponding neighbour pairs, and are given by the 
simple Cluster model as, 


/a-a 

oc 

(1-z) 2 

CQ 

} 

•4= 

oc 

2x(l — x) 

1 B-B 

oc 

Jl 

X 


where x represents impurity concentration. The ratios of the intensities thus calculated have 
a fairly good agreement with the experiments done on KNixMn^xFa except for the Mn — Mn 
band which was found to be much weaker than what is predicted by the Cluster model. This 
discrepancy arose because of a nearly resonant Raman process for the intermediate states 
associated with the virtual excitation of Ni ions. This renders a considerable intensity to the 
Ni - Ni band and even to the Ni - Mn band, since it depends on the overlap of the Ni excited 
states on neighbouring Mn. It is reasonable to assume that the intensity of the Mn - Mn 
band will be considerable only if one or both of the Mn - Mn ions are nearest neighbours 
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t° a Ni ion. So, essentially it is roughly proportional to the Mn pairs which have Ni as 
nearest neighbour. But on increasing the Mn concentration the intensity of the band does 
not increase, as expected naively, due to a simultaneous decrease of nearest neighbour Ni ions. 


The widths of the two-magnon bands are studied in detail by Lockwood et a/.[30]. Eval- 
uation of the two-magnon scattering intensity involves a summation over states within the 
BZ and hence there is a finite linewidth in the scattering spectrum even in the pure materials. 
In presence of impurities there will be two contributions to the width: one arising from the 
wavevector sum and the other from the width of the intrinsic excitations which presumably 
arises due to impurity-induced damping in the magnon spectrum. 

1.4 Experiments done on Copper- Oxide based materials 

Ever since the discovery of high-T c superconductivity[31], there have been persistent efforts 
among the theoreticians to understand the mechanism responsible for pairing among the elec- 
trons. The widespread belief is that the strong electron correlation in these compounds may 
contain a clue to that. One of the important consequences of this correlation is that the 
undoped parent compounds are antiferromagnets. Raman scattering has proved to be a very 
useful tool to probe the collective magnetic excitations in these compounds. Common to all 
antiferromagnetic cuprates, the Raman studies reveal a well defined peak which is attributed 
to inelastic scattering from two-magnon excitations in the Bi g scattering geometry along with 
a weaker but still significant signal in the A lg geometry at a transferred photon frequency 
close to 3000 cm -1 . 

These antiferromagnetic compounds are well described by modeling the undoped Cu0 2 
layer with the Cu ions arranged antiferromagnetically in 3d 9 configuration with the 3 d holes 
carrying a spin 5 = A more detailed picture yields that the Cu spins are coupled to each 
other by superexchange via the Oxygen atoms. 
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In these systems the interaction between the spins may be expressed in terms of the 
Heisenberg Hamiltonian, 


H=j£Si.Sj (1.2) 

( 0 ) 

where S; represents the spin vector at the site i and J is the exchange constant. Conventional 
spin-wave theory yields, 

u,, = ZJS/l - 7, 2 

for the spin-wave spectrum where Z = 4 (for 2D) and S = \. Neutron scattering experi- 
ments on these systems have confirmed the presence of 2D antiferromagnetic order. But due 
to energy limitations in the neutron scattering experiments it was not possible to measure 
the spin-wave dispersion with any degree of accuracy and hence the value of J could not be 

determined. 

So to ascertain the value of J and to probe the spin-dynamics in detail, inelastic light 
scattering studies were performed in these systems and here I am going to provide a brief re- 
view of the experiments done on the insulating and doped phase of the Cuprates. In contrast 
to neutron scattering, these experiments show striking deviation from classical spin-wave be- 
haviour. 

The first experimental observation for the two-magnon Raman spectra came about in 
1987(32, 33], when scattering measurements were performed on antiferromagnetic La 2 Cu0 4 . 
The lineshape is shown in Fig. (1.2). The symmetry consideration plays a crucial role in 
interpreting the polarization dependence of Raman scattering. From the spectral lineshape it 
can be noticed that, 

(a) The two-magnon peak occurs at an energy close to 3 J, to be precise at 2.7J[34] which 
is predicted by a naive bond counting argument which goes as follows. The process of excit- 
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ation of two magnons in a Neel state involves breaking of 6 nearest neighbour bonds which 
results in an energy cost of 3 J as opposed to 4 J, if the spins flipped were not neighbour pairs. 

(b) The anomalously large linewidth of the profile which is fou n d to be approximately 38% of 
the frequency at the spectral maximum, whereas the theoretically predicted value is only 11%. 

(c) There is an appearance of signal in the Ai g and B 2g scattering geometries which are for- 
bidden in the Neel ground state. 

(d) The scattering intensity persists well beyond 4 J whereas the joint magnon density of 
states vanishes at energies higher than 4 J. 

(e) From the experiments the value of J was found to be (1030 ± 50)cm -1 . 


Initially [35] the deviation from the theoretical prediction was suggested to have extrinsic 
origin, such as impurities or non-magnetic scattering (e.g. from phonons). Almost immedi- 
ately the problem was looked at by Singh and co-workers [34]. They took into account the 
effect of quantum spin fluctuation which leads to a substantial enhancement of the linewidth. 
Also the signals in the classically forbidden scattering geometries e.g. A lg and B 2g can be ar- 
gued to originate from diagonal nearest neighbour (NNN) or next nearest neighbour (NNNN) 
spin pair excitations which are Raman active solely by virtue of quantum fluctuations. To 
account for such a NNN term (or NNNN term) in the Raman operator, (see Eq. 1.1) S has to 
be changed to S' where S' connects to NNN (or NNNN) spins. Such a replacement makes no 
contribution to the Bi g geometry but accounts for the spectra in the Ai g and B 2g geometries. 

Also the temperature dependence of the Raman spectra is investigated[35] and it is 
observed that the renormalisation of the overall Raman intensity with temperature is quite 
small compared to that of long wavelength magnons, because of the dominance of short 
wavelength magnons in the two-magnon scattering process. Typically the spectral weight 
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Figure 1.2: Two-magnon Raman spectra for La 2 Cu 04 in the (a) Bi g and (b) Ai g and B 2 g 
scattering geometry. In (a) the peak frequency cc p (= 3170cm -1 ), is found to be slightly smaller 
than the first moment or the ’central frequency’ p\. The second cumulant M-i measures the 
width of the peak. In (b) cr, a' are unit vectors connecting to the nearest neighbour (NN) 
and diagonal nearest neighbour (NNN) in a lattice. (From Ref. [34]) 


shifts towards lower energies by less than 5% between T = 0 and T — T/v(14]. 
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Figure 1.3: Two-magnon Raman scattering intensity shown for YBa 2 Cu306 as a function of 
transferred photon frequency. (From Ref. [44]) 

Similar features are obtained for YBa 2 Cu306[35, 37], The two-magnon peak occurs at 
somewhat lower in frequency and hence yielding a value of J ~ 950cm -1 . The linewidth is 
even greater than La 2 CuO-i (Fig. (1.3)) 

Finally, the Raman scattering cross section was investigated in doped cuprates. Sugai 
d a/. [38] and Sugai[39] obtained data on two-magnon Raman scattering in (Lai_ x Sr x ) 2 CuO^ 
for various doping concentration ranging from 0.0 to 0.06 (Fig. (1.4)). It is observed that 
the position of the peak shifts towards lower frequencies and concomitantly the spectrum be- 
comes broader. The peak structure almost disappears at x ~ 0.06. The scattering spectrum 
in YBa 2 Cu 30 < 5 +x is more robust to increase of oxygen concentration. The two-magnon peak 
survives upto x ~ 0.4 before it becomes almost flat. The rapid suppression of the spectral 
weight is still an open issue. It could probably be related to the loss of intermediate charge 
transfer states of the two-magnon scattering process. This suggestion is consistent with the 
optical conductivity expcriments[40] which demonstrates the transfer of spectral weight from 
above the charge transfer gap to the Drude peak. It was attempted by Lyons and Fleury[41] 
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to explain the effect of carriers on tire spin-excitation spectrum using simple ideas depend- 
ing upon whether the carriers reside predominantly on the Cu 2+ or O 2- site. In the former 
case, Cu 3+ , having no spin, appears as a vacancy (non-magnetic impurity with, of course, no 
potential barrier locally) in the lattice. In the latter case, the O - acquires a spin and thus 
can interact with nearest neighbour Cu 2+ ions. Regardless of the sign of this interaction, 
Cu — Cu interaction turns ferromagnetic. This sign reversal in J is a local phenomenon and 
may have a more disruptive effect upon antiferromagnetism than a simple absence of spin in 
the Cu 3+ case. But in either of the situations, disorder dilution of the lattice occurs. These 
ideas do not treat the holes as mobile vacancies. Till date, there is no firm theoretical pre- 
diction about the effect of carriers on the Raman spectra. 



Figure 1.4: Two-magnon Raman spectra of YBa 2 Cu 3 0 6 4 -x for three values of x: (a) 0.0; (b) 
~ 0-7; (c)^ 0,9. Incident wavelength is 4880A 0 . (From Ref. [37]) 


The Floury- Loudon theory works well for the incoming and outgoing photon frequen- 
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to be considerably smaller than the energy gap between the conduction and valence bands 
in copper oxide insulators. But in reality, Raman scattering experiments done on these 
cuprates, the frequencies are close to the energy gap value (both lying between 2 - 3 eV). A 
comprehensive study to clarify quantitatively the resonance behaviour was initiated by various 
groups[42, 43]. Blumberg et al. reported resonant Raman lineshape for single CuC >2 layer 
Si' 2 Cu 02 Cl 2 and bilayered antiferromagnetic systems like YBa 2 Cu 30 e+ x with a frequency 
range from 1.65 eV to 3.05 eV. Their data shows a composite structure of the lineshape and 
a strong non-monotonic dependence of the scattering intensity as a function of the incoming 
photon frequency. The results were analyzed using the idea of ’Triple Resonance’ due to 
Chubukov and Frenkel[44]. 

1.5 Recent Theoretical Developments 

In the early nineties, Shastry and Shraiman[5] have developed a systematic formalism for un- 
derstanding Raman scattering from Mott insulators starting from an electronic Hamiltonian. 
In the absence of any resonance enhancements, they derived the Fleury-Loudon Hamiltonian. 
However, they find that, when one considers the scattering with the photon frequency close 
to the band gap in cuprates, in principle a large number of parameters are needed to describe 
the scattering process. One can organize the various terms by irreducible representations 
of the lattice symmetry group and by powers of t/(U — cjj), where t represents the hopping 
matrix element, U, the on-site Coulomb term and o\- being the frequency of the laser beam 
used. Within this theoretical framework, the intensities in various geometries are, in general 
not fixed and it may so happen that the signal in Ai g (or B 2g ) geometry is stronger than that 
in Bi g geometry provided u;,- is appropriately tuned. 

An interesting new term was observed by Shastry and Shraiman which is a scalar triple 
product of three spins at different lattice sites, which represents scattering from chiral fluc- 
tuations. They have also argued that the scattering of circularly polarized light can not be 
deduced from linearly polarized light, and by comparing various scattering intensities, in- 
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rluding the circularly polarized ones, one can disentangle the contributions from this chiral 
term. 

To focus on the main aspects of the calculations due to Shastry and Shraiman, one 
can write down the Hubbard Hamiltonian which includes the interaction of light with the 
electronic system as a perturbation, 

H = Ho - t Y2( c l c J^ ie/hcfA - dl + h.c.) 

<d> 

where Ho is the unperturbed part, 

Ho = U £ n, t n i4 . 

t 

and A" is the n-th component of the time dependent vector potential at wavevector q. Taking 
a Fourier transform allows one to write, 

H = Ho - e£jf At, + -Y, 

9 9l .92 

The component of current J* and the inverse mass tensor are defined as, 

ja Y' t 

J q 2-^ flL C fc+a a C k-%,cr 

k,a 

and, 

a 8 SP ^ t 

T ' _ S dk a dk 0 

The Raman cross-section can now be obtained using Fermi’s Golden rule where the initial 
and the final states are direct product of one-photon states and appropriate matter states (i.e. 
states of the Hubbard Model) with energies Ei = u + w* and Ej = e f +u> f . So the inelastic 
scattering rate for transitions can be expressed as, 

r = E|(/l*K)|h(£/-a) 

i,f 

where M is the Raman scattering operator comprising of two terms - a resonant (Mr) and a 
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on-resonant term (M,\). Hence, 

</ivio = + M°/\i) 

'f and (ey)’ are the polarization of the incoming and outgoing photons. M'jf is the reson- 
mt term which is relevant for discussion of the insulating phase of cuprates, whereas the 
ion-resonant term becomes important for the doped case, since the inverse mass tensor 
r 9 °i +92 ^ as non ~zero matrix elements between the states in the lower band, i.e. with low energy 
transfer. 


The resonant term can be written as a combination of two terms — a direct and and an 
exchange term shown in Fig. (1.5), 


{/|A/ fi |0 = 53(0; cr / |j M (r)|l;a)(l;a|;;(r / )|0;cr,) 


pi p* 


+ 


P x pi 


[ U — U){ U + U)f 


M ^ can be decomposed into various scattering geometries pertaining to the D 4 h tetragonal 
point group. For convenience of the future discussion three of the most commonly used 
scattering geometries are shown in Fig. (1.6)(46|. It was shown by Shastry and Shraiman 
that Mr can be expressed, to lowest order in ^ , as the Fleury-Loudon Hamiltonian. 

As mentioned earlier, the expansion in the Shastry and Shraiman was actually in powers 
of — — (when u>, -C l then the expansion can be thought to be in power of jj as mentioned 
at the end of last paragraph) and for the light scattering experiments on high-T c cuprates u is 
roughly of the order U. Then one needs to question the validity of Fleury-Loudon Hamilto- 
nian to calculate the scattering cross-section. 
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Figure 1.5: Diagrammatic representation of the Direct (top) and exchange (below) terms 
appearing in the resonant Matrix element Mr 3 (see the discussion in the text) 


Chubukov and Frenkel[44] studied the case when in; or (w/) ce U (resonant two-magnon 


scattering) and looked at the variation of the overall scale of the coupling constant A (see 
Ec,. 1.1) with the incident photon frequency in. The problem was tackled diagrammatically 
and only the diagrams (Fig. (1.7)) which contribute to the intraband processes are retained. 


This idea has previously been applied for calculating the Raman intensity from two-phonon 
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Figure 1.6: Commonly used scattering geometries. 

process. [47] From the explicit frequency integral of the diagrams, it can be noticed that these 
diagrams have more number of denominators than the ones relevant for the Fleury- Loudon 
case. There are in fact three of them, and all three vanish simultaneously for the photon and 
magnon frequencies satisfying some specific energy-momentum relation. This phenomenon is 
called as ’Triple Resonance’ by the authors. Also they find that the ’Triple Resonance’ yields 
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the strongest divergence of the Raman scattering cross-section and alongwith incorporating 
the spin density wave (SDW) dispersion relation for the carriers, the unusual experimental 
features of the two-magnon lineshape and the two-magnon peak intensity dependence on the 
incoming photon frequency can be explained. 

The success of the theory relates to the experiments recently done on the Cu-Oxide based 
materials. The two-peak structure obtained in the two-magnon profile in Pr 2 Cu 04 [ 4 S], a 
stronger one at ^ ~ 3 J and another one at ~ 4J is in perfect accordance with the theory. 
The theory predicts scattering in both Bi g and Ai g geometry which is an experimental reality, 
as seen before. Furthermore, the calculation supports a single maximum, while the diagrams 
relevant to the Fleury-Loudon case might have given two peaks, one each at incoming and 
outgoing photon frequencies, commonly termed as incoming and outgoing resonances. [49]. 




(b) 


Figure 1.7: Diagrams which are important for the resonant processes. All the three denom- 
inators in the diagram (a) can vanish simultaneously, which is known as ’Triple Resonance’. 
Additional relevant diagram (which is less dominant than (a)) is shown in (b). Here the 
fermions belonging to a valence (conduction) band is represented by a dashed (solid) line. 
The emitted magnons are denoted by the solid wavy lines whereas the dashed wavy lines are 
the incoming and outgoing photons. (From Ref. [44]) 

Recently, Schonfeld et a/. [50] revisited the issue of resonance enhancement of the two- 
magnon intensity in a spin-density-wave (SDW) antiferromagnet. In their diagrammatic 
formulation they have carried out a systematic microscopic derivation for the two-magnon 
scattering intensity explicitly taking into account the resonance structure of the photon-- 
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magnon coupling vertex. Essentially they confirmed the conclusions of the Triple Resonance 
theory due to Chubukov and Frenkel for the resonant scattering process. 

Morr et a/. [51] have studied two-magnon Raman scattering in a two-layer antiferromag- 
net for both Fleury- Loudon and the resonant case. It enables one to obtain an estimate of the 
interplanar exchange constant (</jl). In particular, they observed qualitatively new features 
for the scattering profile in the A lg geometry. It is easy to show that a non-zero J± gives a 
signal in A lg geometry even for the Fleury-Loudon case. In addition to that, they obtained 
a strong resonance in the A ig geometry for the transferred photon frequency ui (= u.', -u s ) 
to be equal to 4(J||Jx) 1 ^ 2 and at resonance, the amplitude of the Ai g signal is larger than 
that of the B lg at the same frequency. The value of Jj. as obtained from their calculation is 
~ 0.1 Jy and the upper limit is ~ 0.25 Jy. Hence J L for the double-layered YBa 2 Cu 3 0 6 can 
be estimated to be 12 MeV < J x < 30 MeV for Jy ~ 120 MeV. This is consistent with the 
NMR data[52] 5MeV < J L < 20 MeV obtained on double-layered Y 2 Ba 4 Cu 7 0i5, but has 
poor agreement with the infrared transmission and reflection measurements [53] where Jj_ is 
0.55 J\\. 

Recently Nori et a/.[54] took into account of magnon damping arising due to the in- 
teraction between the magnons and the lattice vibrations or phonons. They calculated the 
Raman spectrum for a S = 1/2 Heisenberg antiferromagnet by Exact Diagonalization and 
Quantum Monte Carlo techniques for systems ranging from 16 to 144 sites. They propose 
that the phonon-magnon interaction leads to a random fluctuation of the Heisenberg exchange 
J. They assumed that the random couplings SJij can be drawn for each bond independently 
from a Gaussian distribution of width a. A simple calculation to get an estimate of a leads to 
a rather high value ~ 0.4 J. Such an enormous variation of the exchange coupling seems phys- 
ically unreasonable[50]. However, the Raman spectra calculated taking this random model 
show good agreement with the experiments done on the undoped cuprates. The lineshape 
shows a substantial broadening, a profound asymmetry appears in the higher frequency re- 
gion due to the existence of a spectral weight and the intensity persists well beyond 4J. In 
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Fig. (Tb), the lineshapes in the Bi g and Aj g geometries as obtained by Nori are shown 
which show a strong evidence that the Raman lineshape gets considerably modified due to 
interaction with the low-energy phonons which randomly distort the lattice. 



Figure 1.8: Normalized Raman cross-section /(u>)// mox is shown as a function of scattered 
photon frequency c o (in units of J). The plot in solid line shown in (a) corresponds to the 
spectra obtained by Exact DJagonalization with N = 16 and the dashed one is obtained from 
Quantum Monte Carlo studies for N = 144 for a — 0.0 in the Bi g symmetry. The plots shown 
in (b) and (c) corresponds to the scattering intensity in B| g and A lg symmetry corresponding 
to a = 0.3J, 0.4J, 0.5J. (From Ref. [54]) 
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Also Canali and Girvin[55] have attempted a systematic analysis within the spin-wave 
theory to incorporate the contribution coming from the four-magnon process to explain some 
of the anomalous features mentioned in the last paragraph. Although the first three moments 
of the scattering intensity are in good agreement with experiments and the series expansion 
results due to Singh at al. [34] , but the full lineshape could not be successfully explained. 
The integrated four-magnon intensity is merely 2.9% of that coming from two-magnon pro- 
cesses. So, it seems the four-magnon contribution can barely modify the overall two-magnon 
Raman lineshape. Furthermore, they observed that the four-magnon peak frequency is 2.5 
times larger than that of two-magnon. Hence this leads one to conclude that the four-magnon 
scattering can not contribute to the spectral intensity at a frequency close to 4 J. 

Finally, before winding up the discussion on the theoretical progress that has taken place 
in recent times, I present a brief account of the Raman scattering studies performed in the 
t — J model in presence of holes[56, 57]. The calculations were performed using Exact Diag- 
onalization studies on a 4 x 4 lattice in the B lg scattering geometry. The cases arising out of 
no hole, one hole and two holes are discussed. The study for the no-hole case reveals good 
agreement with the series expansion results obtained for the cumulants of the distribution[34]. 
In presence of holes the two-magnon shifts peak downward in energy and the intensity reduces 
with respect to the no-hole case, while the spectral weight appears at higher energies. The 
position of the two-magnon peak for the one hole (static) case as a function of transferred 
photon frequency has got good agreement with what one expects in the Ising limit from the 
flipping of two spins in the vicinity of the hole where the number of broken bonds created by 
the Fleury-Loudon term is only five while away from the hole is six. For the two-hole case 
(where the two holes are neighbours in the ground state in order to minimize the number of 
broken bonds) 3 possibilities arise, viz. both the spins flip away from the holes, one of them 
in the vicinity and both are in the vicinity of the holes. In the Ising limit the peaks occurring 
from such cases should be separated in energy by ~ 0.5J and are in good agreement with 
the numerical work. 
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For a dynamical hole, Gagliano and Bacci[57] find that the spectrum broadens out (see 
big. (1.9)) vvhile the two-magnon peak at ~ 3 J does not continue to be the dominant one 
and more peak structures appear in the lower energy regime. These effects have been experi- 
mentally observed[38] and show the relevance of doping concentration on the spin-fluctuation 
dynamics. 



Figure 1.9: B lg Raman spectrum for J = 1 in a 4 x 4 lattice, obtained from Exact Diag- 
onalization studies, (a) Pure Heisenberg model; (b) static vacancy in the Heisenberg, (c) 
dynamical hole in the t - ./*; (d) dynamical hole in t - J model. (From Ref. [57]) 
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In spite of a great deal of success in explaining the two-magnon Raman lineshape for 
the insulating cuprates, at least three characteristic features do not get completely resolved - 
anomalously large linewidth of the profile, asymmetry towards the higher energy region with 
respect to the two-magnon peak and the lineshape considerably extending beyond 4 J which 
corresponds to the maximum of joint magnon DOS. In my thesis, I have addressed these 
issues and in addition, investigated the effect of disorder on two-magnon Raman scattering 
spectra in arrtiferrornagnets. 

1.6 Brief Outline of the Thesis 

In Chapter 2 we have developed a systematic diagrammatic scheme in order to study two- 
magnon Raman scattering in a pure antiferromagnet. There the fermionic structure of the 
magnon interaction vertex is studied at order l/Af where J\f denotes the orbital degeneracy. 
The full two-magnon propagator is then evaluated by including the interaction between the 
magnons at the ladder sum level and the scattering intensity is obtained by taking the imagin- 
ary part of the propagator. Finally the calculations are extended to the intermediate coupling 
regime in order to make connections with the experiments done on cuprate insulators. We 
also main* an attempt to incorporate zero temperature magnon damping in a phenomenolo- 
gical way which yields a better agreement with the experiments performed on La 2 Cu 04 . 

In Chapter 3 we study the effect of hopping disorder in an antiferromagnet using both 
perturbation theory and exact eigenstate method. Here the lattice distortions arising out of 
quantum and thermal fluctuations are taken into account which essentially renders the ex- 
change couplings Jij. and hence the hopping integral tij, which depend upon the instantaneous 
positions of the ions, random. It was shown that while the low-energy long wavelength, low 
energy magnon modes are only weakly affected, the high energy modes are strongly renor- 
malized because of a cooperative effect arising from local correlations in hopping disorder. 
We also discussed the implication of this highly asymmetric magnon energy renormalisation 
for two-magnon Raman scattering in cuprate antiferromagnets. 



1.6 Brief Outline of the Thesis 


30 


Chapter 4 deals with obtaining the Raman scattering intensity numerically in presence 
of hopping disorder. A consistent scheme of mode assignments is prescribed which estab- 
lishes a one-to-one correspondence between the pure and disordered system. It was shown 
that for a substantially weak disorder, several features like the linewidth to peak position 
ratio, asymmetry with respect to the two-magnon peak and the spectral weight persisting all 
the way upto 6 J agree with the experimental data. 

In Chapter 5 we have calculated the Raman scattering cross-section in presence of on- 
site disorder. It is found that the magnon modes are damped and this damping is obtained 
by considering the vertex corrcctions[58]. Even though it is observed that this damping is 
proportional to fourth power of the disorder strength, still for the high energy modes which 
are particularly relevant for the Raman studies, the damping goes as the DOS which shows a 
divergence at the upper band edge (2 J) and hence is a significant quantity. We have utilized 
this damping to calculate the Raman lineshape for different disorder strengths. It is seen that 
there occurs a broadening of lineshape with a noticeable change towards the high frequency 
regime, whereas the correction in the low energy region is negligible. 

In Chapter 6 we develop a formalism to obtain the two-magnon Raman scattering 
intensity in real space. This involves calculation of the two-magnon propagator in terms of 
the eigensolutions of the zeroth order particle-hole propagator [y 0 ]. Finally the imaginary 
part of the two-magnon propagator is taken to obtain the scattering intensity which can be 
computed for various kinds of disorder such as hopping disorder, potential disorder, vacancies 
etc. 
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Chapter 2 

Two-Magnon Raman Scattering for a 
Pure Antiferro magnet 

2.1 Introduction 

As mentioned in the last chapter extensive studies have been carried out in antiferromag- 
nets such as RbMnF 3 , MnF 2 , K 2 NiF 4 ,[l-5] and a good theoretical understanding exists for 
these systems with spin 5 > 1, within the interacting-magnon theory for the Quantum Heis- 
enberg Antiferromagnet (QHAF)[6]. Recently there has been renewed interest since stud- 
ies in La 2 OuO,i (a S = 1/2 system) have indicated evidence of substantial quantum spin 
fluctuations, [7-9] and the observed two-magnon Raman scattering linewidth is much too broad 
to be explained within the classical theory.[4] Moreover, scattering intensities are observed in 
geometries which, in the classical approximation, should yield no scattering. Interest in two- 
magnon Raman scattering continues also because upon doping with holes, the two-magnon 
peak broadens out and relaxes into the Raman continuum. In YBa 2 Cu306+ x , for example, 
the intensity of the two-magnon feature decreases rapidly with increasing x, and virtually 
disappears at x = 0.5. [10] 

While theoretically two-magnon Raman scattering has been studied extensively within 
the QH AF, there have been so far only limited attempts to systematically study this within the 
Mott- Hubbard model. Recently a strong-coupling expansion in powers of t/U , the hopping 
term, has been carried out, and the Fleury- Loudon Hamiltonian has been obtained at the 
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second order level. [11] In this treatment, the hopping term is extended to include the gauge 
teim t c ■ arising horn the external transverse electro-magnetic field, and the transition 
matiix elements which determine the Raman scattering are obtained within a strong-coupling 
expansion. 1 his has also been described in a somewhat detailed manner in the last chapter. 
However it appears that so far no systematic, formally weak-coupling expansion has been 
carried out which can continuously interpolate between the small -U and larg e-U limits. 

Significant!}, in the small- and intermediate-?/ regimes of the Mott-Hubbard antiferro- 
magnet , when the interaction term U is less than or comparable to the free-electron band- 
idth B = 2Zt where Z = 2D, D being the dimensionality, the one-magnon density of states 


w 


(DOS) exhibits a peak structure which broadens and shifts to lower energies with decreas- 
ing / .[12, 13] In contrast, in the large-// limit, the one-magnon DOS diverges at the upper 
band edge due to zone-boundary magnons at energy DJ, where J = 4t 2 /U is the exchange 
energy. This peak structure resembles a broadening due to magnon-interaction, and must 
be incorporated in any realistic theory of two-magnon Raman scattering when applied to 
the small- and intermediate- // regimes. Since it is believed that La 2 Cu 04 does actually fall 
in the intermediate-// regime, a formally weak-coupling expansion scheme has a potential 
application as well. 

In this chapter we describe a perturbation-theoretic, diagrammatic scheme for systemat- 
ically obtaining the two-magnon Raman scattering in an antiferromagnet. Since within the 
Hubbard model itself there is no formal expansion parameter, we make use of the generalized 
Hubbard model with A" orbitals per site, [14] and develop a systematic expansion in powers of 
l /A/ This perturbation-theoretic scheme, which preserves the spin rotational symmetry, and 
hence tin* ( JoldsUme mode, order by order in the perturbation theory, has been used earlier to 


systematically obtain quantum corrections to sublattice magnetization, spin-wave energy, per- 
pendicular susceptibility and ground-state energy of a Hubbard antiferromagnet. [14] When 
this formally weak-coupling perturbative approach is carried to the large-// (strong-coupling) 
limit, which is analytically the simplest, one recovers identical results as obtained within the 
1 /./-expansion for the spin S Quantum Heisenberg model. We shall see that this holds for 
the two-magnon Raman scattering effect as well. 
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In spin-paii excitations by light, the two magnons are created on nearest-neighbour (NN) 
sites, and theiefore stiong interaction effects between magnons are important in a quantitative 
analysis. Considering a two-dimensional 5=1/2 system for concreteness, two spins flipped 
involve an excitation energy 4 J when far apart but only energy 3 J when on NN sites. This 
SU SS (>!:! ^ S that the interaction energy between magnons at the lowest order level is — J when 
on NN sites, which shifts the peak in the two-magnon Raman scattering intensity to energy 
~ 3.7. Indeed, we find within our perturbation-theoretic approach, that to the lowest order 
in 1/.V , the interaction energy between magnons on NN sites is instantaneous and precisely 
— ,/ in the large-// limit. 

In section 2.2 the two-magnon propagator and magnon-magnon interaction are con- 
sidered at the lowest-order level in an expansion in powers of f 2 //7 2 . In the next section we 
have extended the calculations to the intermediate-// regime by incorporating higher order 
terms in the magnon propagator and interaction terms. Finally in section 2.4 we make a 
simple estimate for zero temperature magnon damping and incorporate it to the evaluate 
the Raman scattering intensity in the B lg symmetry and compare with the observed Raman 
spectrum in La^CuO,!. 

2.2 Two-magnon Propagator in the large - U limit 

The Floury- Loudon Hamiltonian which represents the interaction of spin pairs with photon 
pairs, as appears in Eq. (1.1), is written as, 


Hr = A £(E in< J)(E sc i)S(r).S(r + $), (2.1) 

rj 

where the various symbols have been described in chapter 1. The sum over r ensures that 
the spin-pair excitation has zero total wavevector as required by momentum conservation 
since the photons involved have essentially zero wavevector. To calculate the two-magnon 
spectrum, we need to extract parts of the spin operators which combine to create a pair of 
magnons. We focus therefore on the transverse part of the spin pair operator, 
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p } = ES(r).S(rtJ). 

r 

AvS mentioned in the last chapter, the two-magnon Raman scattering intensity is related to 
the con elation function of the spin-pair operator Pj, which can be obtained as the imaginary 
part of the corresponding propagator. We consider the following time-ordered two-magnon 

Green's function, 

0 U ,m = -i / Er 5-(r,i)S + (r + j,i) (2.2) 

£ r ' S + (r',i')5-(r' + i',<')]|V’c> 

As mentioned earlier, in this chapter vve develop a systematic perturbation expansion for 
this two-magnon propagator within the Hubbard model. We show the fermionic structure of 
the magnon-interaction vertex to 0{l/A r ), and evaluate the resulting two-magnon propagator 
in the larg e-f r limit in order to make contact with the known results within the Heisenberg 
model formalism. Extension to higher orders in (1/A/”) and to the intermediate-// regime can 
tlu'ii be carried out. 

At this point we will take a short break to discuss the spin-wave excitations in an 
antiferromagnet[I2]. Spin waves are the lowest lying collective excitations that represent 
transverse spin fluctuations about the antiferromagnetic ground state. The spin-wave modes 
show up in the form of a pole in the transverse spin susceptibility evaluated in the AF ground 
state. The transverse spin susceptibility is defined in terms of the spin raising and lowering 
operators as. 


v - + (r.i;r',0 = i9(i-t')(*AF|[S-(>-,<).S + (r',i')]|^ f .). (2.3) 

Within random phase approximation (RPA) or at the ladder sum level for the above propag- 

ator one obtains. 
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The zeioth-order term in the susceptibility, x°, is evaluated with respect to the AF ground 
state and is given b} the following relation in terms of the one-particle Green’s functions, 

\ 0 ' r, r'\ ft) = i f+“ ^G t {r,r , ] Sl')G i '(r',r-,n' - to). (2.5) 

The full RI A propagator \ + corresponding to the pure case has been discussed elaborately 
earlier in the large- ('[12] and intermediate-£/[12] regime. 

Now we consider the non-interacting limit of the two-magnon propagator, G°- . (ft), 

5, S f 

decompose it into the irreducible representations of the lattice, and then develop the perturb- 
ative expansion in powers of the magnon-interaction. In the non-interacting limit we have 
a magnon pair propagating from NN sites r, r + 6 to r', r' + S', and in terms of the magnon 
propagators, \ + and \ H , we obtain for the two-magnon propagator in the non-interacting 
limit: 


G uW 


i J cite ,n(i_t/) £ x + (r, t; r', t')x + ~( r + S , t ; r' + S', t') (2.6) 


r, r' 


Here we note that since S,S' are vectors joining nearest neighbours, if r and r' belong 
to sublattices ft and u respectively (fx,u = A/B), then r + 6 and r' + S' belong to opposite 
sublattices /7 and u respectively, and therefore we obtain, 




/#EE[x- + W.ft.)Ux + -(-o.n-o,)]^. 


,Q.(S-6‘) 


(2.7) 


Q A'*' 


It is convenient to (i) decompose G\ f (ft) in terms of the irreducible representations 

o, S' 

<p n {Q) of the lattice and (ii) consider G° in a 2 x 2 matrix form within the two-sublattice 
basis. This greatly simplifies the summation of the perturbative series in powers of magnon 
interaction, as we will see later. For the n-th irreducible representation, we obtain, 
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c;° n (n) 


f dCl 

1 j [x“ + (Q,ni)]^ [x + "(-Q,n-ni)]_ (2.8) 


1.1V 


Finite-temperature evaluation of [G°(0)] Ml/ can now be carried out in the standard way. 
Replacing the frequency f l by the bosonic Matsubara frequency iCl m and the frequency in- 
tegial ovoi by a contour integral, and taking the contour around each singularity in the 

clockwise direction, we obtain, 






'7 


d£l i 


E<?n (Q) [x~ + (Q^i)\ uu lx + -(-Q,ittm - no 


J fJLV 


liv g 




(2.9) 


Adding the residues from the four poles, taking e 0tQm = 1, and finally replacing iCl m by 
fl, yields the following matrix elements of [G'°(f2)][12], 




BB 



ErtlQ) 

Q 

Z*UQ) 

Q 


'D£ 

'D£ s 
Qq , 


1 + 

1 + 


Q 




Adding these four matrix elements allows for all possibilities regarding positions of r and r', 
and imaginary part of the 1 sum yields the non-interacting limit of the two-magnon Raman 

scattering intensity. 


iKm = X>t«) (^) V-2««)coth(^£) (2.11) 

which is essentially the joint magnon density of states near the upper end of the magnon 
spectrum where* Uq ~ I)J, and peaks at 0 = 2 DJ. Finite temperature corrections are 

negligible at temperatures k B T < J. So the expression reduces to, 
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= MW) (|£) «(n-2n Q ) (2.12) 

Wo now turn to the magnon-interaction vertex within the Hubbard model, and obtain 
its fermionic structure at the lowest order level 0(1 /M). The magnon-interaction vertex is 
evaluated in the large- f r limit, wherein we only retain terms of order t 2 /U 2 and neglect terms 
of order t A /l A and smaller. t 4 /U A ~ J 2 /U 2 is a measureof the next-nearest-neighbour (NNN) 
spin-spin interaction within the equivalent extended-range Heisenberg model, and 0(J 2 /U 2 ) 
terms will be included in the intermediate - U regime. The fermionic structure of the magnon 
interaction vortex, when the interacting magnons are on opposite sublattices, is shown in Fig. 
(2.1 ). There are 8 diagrams which contribute at the 0(1 /Af) level to the magnon interaction 
vertex. These diagrams essentially describe various f-spin and 4-spin particle-hole processes 
involving the nearest, neighbouring sites. The +/ — signs in these diagrams refers to the 
adcanced/retarded nature of the one-particle Green’s function, corresponding to states in 
the upper/lower Hubbard bands. The fermionic propagator signs indicated in the diagrams 
are not exhaustive of all possibilities, but correspond to choices which yield leading-order 
contribution at the t 2 /(’ 2 level. The contributions of diagrams shown in Fig. 2.1 (a) - (h) are 
given in 'fable 2.1, and detailed evaluation of one of the diagrams is given in the Appendix A. 
Fpto order t 2 /( r ~ the net contribution of all diagrams of 0(1 /Af) to the magnon interaction 

vertex is given by. 


- O', «I - fta) = -1DJ 1q _ Q> (2.13) 

implying nearest-neighbour, instantaneous interaction energy -J between magnons. In the 
magnon-interaction vertex the entering and exiting magnon lines involve the same sublattice, 
and in fact the same site in larg e-U limit, so that in real space we can write 


\ (r,r',r + $,r' + S') = V int (r, r, r + 6, r + S) = - J 
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Therefore in the two sublattice basis the magnon-interaction vertex is diagonal: 

['WQ - (?')]„„ = -ZJiQ _ Q'Kv 

where Z = 2D is the coordination number of the lattice. We also make use of the irreducible 
representations o n (Q) of the lattice to decouple the momentum dependence of the interaction 

vertex as, 

v„AQ - Q') = —-J-IQ -Q' = -JEMQ)MQ') (2.14) 

n 

Hen 1 the irreducible representations <p n {Q) are normalized so that Y,q 4%i(Q) — 1 f° r n - 
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Figure 2.1: Fermionic structure of the diagrams contributing to the magnon interaction vertex 
at the order l/Af level and in the strong coupling limit. The interacting magnons are on A 
and B sublattice sites. The interaction vertices lie within the four wavy lines representing the 
Hubbard interaction. The hatched areas represent magnon propagators. 
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Table 2.1: Contributions to the magnon interaction vertex of diagrams shown in Fig. 2.1 (a)- 
(h). for various signs of the pair of fermion propagators which are off-diagonal in sublattice 
basis. The feimion-piopagator signs refer to the advanced/retarded nature of the fermionic 
propagators. 


Diagrams 

Fermion-propagator signs 

Contribution 

Fig. 2. 1(a), (b) 

(++), ( — ) 

(J D/2)7q_q/ 

Fig. 2. 1(a), (b) 

(+-).(-+) 

-2(JD/2)7q_q/ 

Fig. 2. 1(c), (d) 

( — ) 

( J £>/2)7Q-Q' 

Fig. 2.1(c), (d) 

(+-).(-+) 

-(J D/2)^q-q, 

Fig. 2.1(e), (f) 

C++) 

(J D/2)7q_q/ 

Fig. 2. 1(e), (f) 

(d — }, { — h) 

-(J D/2)7<3_q, 

Fig. 2.1(g), (h) 

( + -).(- + ) 

(JD/2)7q_q/ 


This decoupling permits a straightforward summation of the perturbation series for the 
two-magnon propagator [G„(fi)] involving repeated magnon interactions in a ladder sum. 
The orthogonality of the irreducible representations 4> n (Q ) leads to the following result in the 
two sublattice basis. 


[G n m = 


1 + J [G»(!l)] 


(2.15) 


Again we need to add the four matrix elements of [(7 n (0)] so as to allow for all possib- 
ilities regarding positions of r and r'. If we use A(Cl), B(£l) and C'(fl) to denote the matrix 
elements [6'“(f2)]/M, [Gn(^)]ss anci [G° n {£l))AB/BA respectively which are given in Eq. (2.10), 
then the sum of matrix elements of [G n (fi)] is obtained as below. The imaginary part then 
yields R n (il), the two-magnon Raman scattering intensity for a Hubbard antiferromagnet, in 
the strong coupling limit, at the ladder-sum level. 


( 0 \] — iA — — ~ 9^1 

l_\ LmWW- (i + JA)[\ + JB) - J 2 C 2 


(2.16) 


Within the Heisenberg formalism, using a perturbative approach in which magnon- 
magnon interaction was taken into account in a ladder approximation, this result has been 
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obtained earlier by Solyom.[17] In this work four vertex functions were introduced correspond- 
ing to the four possibilities arising from the incoming and outgoing magnon lines belonging to 
the two sublattices. These vertex functions were obtained, within the ladder approximation, 
in the form of coupled equations, which were then solved and added to yield 
On the other hand, by expressing the two-magnon propagator in the two-sublattice basis as 
2 x 2 matrix, we are able to perform the ladder sum in a very straightforward manner, and 
summing the matrix elements at the end yields the same result. Of course our starting point 
is the Hubbard model, and we have treated the Hubbard interaction term perturbatively 
within a formally weak-coupling expansion scheme, and carrying this expansion to the strong 
coupling limit we have recovered the Solyom result. Our approach can therefore be extended 
to the intermediate!/ and small-H regimes as well in principle. 

Finally we consider the Ising limit of the expression for the two-magnon propagator given 
in Eq. (2.16), in order to clearly see the role of the nearest-neighbour magnon interaction. 
In the Ising limit when the spins can only be flipped completely, the spin-wave (spin-flip) 
energy is ZDJ, jq = 0 and therefore A = C = 0 and B = £ Q ~ n+2 £>j) = 

(3-^ - a^m)- In this limit reduces to the following equation which shows 

the pole at ft = 2 DJ - J corresponding to the two-magnon spin-flip excitations created on 
nearest-neighbour sites: 

£[Gn(ft)Ullsing = Q _ 2 £)J + J 

fXU 

Several extensions of this analysis of two-magnon Raman scattering in the Hubbard 
model, carried out here in the large-17 limit, are possible, some of which are discussed 
below. Higher-order corrections in the inverse- degeneracy expansion in powers of IJM may 
be included in the magnon- interaction vertex. Magnon damping which may be systematically 
incorporated in the theory in this manner is of particular interest because, after all, one 
magnon states are not exact eigenstates of the antiferromagnet, and hence damping of magnon 
modes is expected, particularly of the short-wavelength, high-energy modes. Also the one 
magnon lines themselves may be renormalized to include the magnon-interaction effects. 
At the 0(1 /AO level the spin-wave energies are modified by the momentum-independent 
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multiplicative factor Z c ( = 1.16 in 2D e.g.), and the spin-wave amplitude is reduced by the 
same factor as the sublattice magnetization relative to the HF value. Also of interest is the 
extension to the in termed i ate- V regime which is discussed below. 


2.3 Extension of the calculation to Intermediate- [/ re- 
gime 


As discussed in section II, in the intermediate - U regime the Hubbard model maps to an 
extended-range Heisenberg model and the NNN spin couplings of 0 (t 4 /U 3 ) can be obtained 
within the RPA analysis. [12] Alternatively the magnon-interaction vertices of Fig. (2.1) can be 
evaluated upto the next order, 0 (t 4 /U 4 ). Thus with the one-magnon propagator \ -+ (fl) and 
the magnon-interaction vertex both evaluated upto 0 (t 4 /U 4 ), one can systematically study 
the two-magnon Raman scattering in the intermediate-C/ regime. As mentioned earlier this 
regime is of interest because the modification of spin-wave spectrum results in a significantly 
different form of the rnagnon density of states. Instead of the divergence at energy 2 J in 
DOS due to zone-boundary magnons, one instead gets a peak structure with a progressively 
increasing broadening as 1 1 ft decreases. This peak structure resembles the effect of rnagnon 
damping and therefore any quantitative treatment of two-magnon Raman scattering in the 
intermediate-!/ regime must take this intrinsic broadening into account. Here we consider a 
simple extension to the intermediate- 17 regime. 

There are essentially two modifications that are required in our earlier analysis. Both 
the NN magnon-interaction energy and the spin-wave propagator acquire an order t 2 /C 2 
correction. Using results of reference [12] where spin-wave properties were obtained via a 
systematic expansion in powers of t 2 /U 2 , we obtain after a zero-temperature evaluation of 
the two-magnon propagator the following matrix elements of [G° (D)], the zeroth order two- 
magnon propagator, given in Eq. (2.10), 


a(q) = 5>v;($) 

Q 





n i / fif 
_a 2J + 2 \2j) 


ft 2 - 4 Qq 



2 3 Extension of the calculation to Intermediate-?? regir 


B((l) = X>V’(Q) 


C(f)) = 

Q V 1 QJ L 


f %,\ 2 


+ Q iL + lfjI) 2 1 ^ „_ 

2 J 2 \2JJ J Q 2 - 49J q 


f? 2 - 


(2. IS) 


whoie ni 1 It /A is the magnetization at the HF level, 2A = mU is the Hubbard gap, 
and upto order t 2 j A 2 we have the following results in the intermediate-?/ regime: 


1 1q ) — (^ + 3 cos Q x cos Q y — 97q^ 

t 2 / 3 \ 

1 “ A 2 V + "> C ° S C ° S + J 


(2.19) 


A simple estimate for the NN magnon-interaction energy is now made as follows. We have 
seen that in the large-?’ limit, this interaction energy is nothing but the NN spin interaction 
energy. Since this result is consistent with the simple bond counting argument, we expect this 
to hold in the intermediate-?? regime as well. Therefore we take the NN magnon interaction 
energy to be simply the modified NN spin interaction energy. An elementary quantitative 
analysis which relies on a comparison of spin-wave energy forms for the intermediate-?? 
Hubbard model (evaluated upto order t 2 j A 2 ) and the NNN Heisenberg model is given in 
Appendix B and yields for the NN spin interaction energy: Jnn — J (l — 

The two- magnon Raman scattering intensity for the two-dimensional system is now eval- 
uated in the intermediate-?/ regime from G„(f?) = (?°(fi)/[l + JjVivG^f?)] and is shown in 
Fig. (2.2) for the case A /t = 3.5 which corresponds to U ~ B for the two-dimensional case 
where H denotes the bandwidth. We have taken the the symmetry factor as, 

&n(Q) = (COS Q x COS Qy), 


appropriate to the B lg scattering geometry. We make a comparison with the large-?,' (A ft = 
100 ) limit result and it is dear that while there is a 25% increase in the linewidth in going to 
the intermediate ?/, this increase is not sufficient to account for the zero-temperature linewidth 
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Figure '2.2: Calculated two-magnon Raman scattering intensity in the Bi g symmetry in the 
large-C limit (dashed) and intermediate-!/ limit (solid). 

observed in La-jCuO.i. Therefore while this contribution must be taken into account in any 
quantitative analysis, the large broadening seen in La 2 Cu 04 must find an explanation else- 
where. In the following we consider an important source of broadening which is particularly 
relevant, for La^CnO-i which is a low-spin and low-dimensional system. 

2.4 Zero- temperature magnon damping 

It is well established now that quantum spin fluctuations play an important role in L^CuCXt 
and other cuprate antiferromagnets which are spin 1/2 and almost two-dimensional systems. 
There is a substantial reduction in the zero-temperature sublattice magnetization relative to 
the UK (Neel) value. Thus while long-range AF order is maintained in the 2D antiferromag- 
nut, there is substantial amount of spin disorder. A consequence of this fluctuation-induced 
disorder is that spin- wave states are not exact eigenstates of the antiferromagnetic ground 
state, so that spin-wave or magnon damping must necessarily be present. And since two- 
magmm Raman scattering is related to the imaginary part of the two-magnon propagator. 
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magnon damping should play an important role in the observed broadening. 

\\ hile finite-temperature magnon damping in the Heisenberg antiferromagnet has been 
studied in detail, [lb-20] zero-temperature magnon damping is not easy to obtain in the 
anti ferromagnetic insulator. The second-order magnon interaction process does not yield 
any damping at zeio temperature owing to phase-space restriction. This process involves 
a. magnon dot ay ing into three magnons which subsequently reunite. From simple energy- 
moment um t onsei \ation considerations it is seen that for light-like linear dispersion this 
process has a vanishingly small phase space. This argument in fact holds for any convex 
energy-momentum dispersion, as for magnons. 

We discuss below a simple estimate for zero-temperature magnon damping which is based 
on the following disorder analogy. Because of the substantial transverse spin fluctuations 
present in the S = 1/2, D = 2 anti ferromagnetic insulator, magnons propagating in the 
system will see a fairly disordered magnetic lattice in which the spins, while maintaining 
overall AF long-range order, are slowly fluctuating in space and time. These fluctuations are 
slow on the time scale of (2J) -1 and therefore zone-boundary, high-energy magnons of energy 
2./, which are of interest, will see an essentially static disorder. This disorder-analogy picture 
suggests that magnon damping will be proportional to the magnon density of states, which 
results in a qualitatively correct picture — substantial damping for zone-boundary magnons 
and very small damping for low-energy, long-wavelength magnons. Since we interested only 
in the zone-boundary magnons, in view of the divergent magnon density of states at the 
upper edge at 2,1, a self-consistent evaluation of damping is therefore essential. We obtain 
the self-consistent magnon damping for zone-boundary magnons with energy Q, ~ 2J as given 
below: 

T = 7 [x" + (Q,H)] 

<3 

= 7 Tf^ ( fi_n< J )* + p (2 ' 20) 


where 7 is an effective disorder strength which measures the degree of spin disorder due to 
quantum transverse* fluctuations, which we now proceed to estimate. An appropriate measure 
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of imm spin fluctuations is the equal-time correlation function {S~S+), and therefore 
the simplest estimate for the scattering vertex is J(S~S+). And since the disorder strength 
represents two scattering processes, we obtain 7 ~ J 2 (S~S + ) 2 . Now for the D = 2. 5=1/2 
antiferromagnet, we know (S~ S+) = Eq(1 I - l)/2 « 0.2, which is precisely the 
quantum-lluc tuation reduction in (6-) or the sublattice magnetization. 



Q./J 

Figure 2.2: Calculated two- magnon Raman scattering intensity in the Bi g symmetry with 
a magnon damping of V/J = 0.15, and few points from the observed Raman scattering in 

Lu’jCuO.i (taken from ref. 8]. 

Substituting this value of 7, the effective disorder strength, we self-consistently solve 
from Kq. (2.20) the effective one-magnon damping for zone-boundary magnons, and obtain 
V/J 0.15. Tims even for zone-boundary magnons the resulting damping is small — only 
about 7% of the energy. We now incorporate this magnon damping arising from quantum spin 
fluctuations into the two- magnon Raman scattering evaluation from Eq. (2.16) in the large-C 
limit. This is done by simply adding the imaginary part 2iT to the two-magnon energy 2 VIq 
in the energy denominators in Eq. (2.10). We take the symmetry factor 4> n {Q) = (cos Q x - 
cos Q,j), appropriate to the Bi g symmetry, for comparison with Raman scattering observations 
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on La-jC’uO.]. I ho result of this evaluation, together with the Raman intensity observed in 
LajC'uO.j taken from reference 8 is shown in Fig. (2.3). We adjust the energy scales so that 
the observed Hainan peak at around 3200 cm -1 occurs at 3 J where our calculated Raman 
scattering intensity peaks. As we have shown earlier, there will be additional broadening if 
the ('valuation is done in the intermediate- 17 regime, appropriate to the cuprates. It is clear 
that while the linewidth compares reasonably well with the observed FWHM of about 1200 
cm -1 , there is a marked asymmetry in the observed Raman line, and also a significantly 
larger scat (('ring persisting at energies above 47. It has been suggested that four-magnon 
excitations an* responsible for these features. 
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Appendix A 

In this Appendix wo discuss evaluation of diagrams shown in Fig. (2.1) (a)-(h) representing 
magnon interaction vertices relevant for the two-magnon propagator. The vertices are shown 
for interaction between rnagnons on opposite sublattices A and B. The evaluation is performed 
in the st rong coupling limit in which case the interaction is actually effective only on nearest- 
neighbour sites. As mentioned earlier the interaction vertices represent spin-|, spin-4 particle- 
hole processes in which energy-momentum is exchanged by the interacting rnagnons. All 
diagrams shown are of order 1/A'’ in an inverse-degeneracy expansion. 

In each diagram there are precisely two fermionic lines connecting an A-sublattice site 
to a B sublattice site, and therefore off-diagonal in the two-sublattice basis. This pair of lines 
already contributes a factor < 2 /{’ 2 , and therefore we need to retain only the leading order 
contributions from all other terms in each diagram. As a representative example, we discuss 
the evaluation of the diagram shown in Fig. 2.1 (g), where the fermion lines constituting 
the loop are off-diagonal in sublattice basis. The amplitudes of all other Green’s functions 
an 1 taken to he 1 and poles at energy - ±U / 2, in order to obtain the leading-order 
contribution. ( ’ontrilmtions of all diagrams shown in Fig. 2.1 (a)-(h) are given in Table 2.1. 
To the magnon interaction vertex the contribution of diagram shown in Fig. 2.1 (g) is given 
by: 

WM-Q'Mi -^Mb) 


«(»n n 


J :p 12 (fc, n) AA G^{k - Q + Q', ft - fix + )aaGI ( k-Q,n - fli )aa 
[ y (k' + Q,n'-n + ( fc/ + Q',n'-n + n 2 ) B B 

J 2t7T 

[ l B22 Y tt')A8G${k" -Q + Q', ft" - fti + ^2 )ba 

J 2 7T . ,, 


777 ^ - {~U/2) - iff/ - UI2 + in 
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ft — U /2 + ir/ J 

( e J£_ -Q + Q^_ \ 
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^ (- 1 ) 
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Appendix B 


We discuss below an estimation of the NN spin interaction energy in the intermediated 
regime. I he idea is to compare the spin-wave energy expressions obtained from (i) the 
spin- 1/2 Heisenberg model with NNN spin interaction, and (ii) the intermediated Hubbard 
model, both for the two-dimensional, square-lattice case. If J and J' refer to the NN and NNN 

Heisenberg antiferromagnetic spin interactions, respectively, then for the spin-wave energy we 

have: 


« 2./[(l - 7 2 ) - (2J'/J)(l - 7 ^)] 1 /2 (2.22) 

when' y Q » cos Q r cos Q v% and terms of order (J'/J) 2 have been neglected. We compare 
this expression now with the result for the Hubbard model obtained at the RPA level by an 
expansion[ 1‘2] in powers of t 1 j A 2 , where 2A = mU is the Hubbard gap. 

= 2(4/V«0[(l - 7 1) ~ (t 2 / A 2 )(6 + 3j'q - 97J)] 1 / 2 (2.23) 

which, after separating the 7 q and 7 q pieces, can be rewritten as: 

Slg a 2{4I-/U)(l - 9l7A 2 ) 1/2 ((l - 7 2 ) + (3<7A 2 )(1 - (2.24) 

Now comparing these two spin-wave energy expressions, we obtain for the NN spin 

coupling: 


J 




and a ft rronuujnt tie NNN spin coupling: 


(2.25) 
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Appendix C 


In this appendix we present, a list of the irreducible representations for the square lattice 
antiforroinagnet which are relevant for our purpose, 


MQ) = 


v/2 


(cos Q x COS Q y ) 


(2.27) 


MQ) = 



cos Q x d~ cos Q y) 


(2.28) 


<h[Q) = 


_ 1 _ 

72 


(sin Q x 


- sinQy) 


(2.29) 


MQ) 


V 2 


(sin Q x + sin Q y ) 


(2.30) 

(2.31) 


Among these, the first two are appropriate to the Bi g and A ig scattering geometry. Also the 
orthonormality property of these irreducible representations can be written as, 


= 1 for a11 n 

Q 
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Chapter 3 

Hopping Disorder and Asymmetric 
Magnon-Energy Renormalization in an 
Antiferromagnet 

3.1 Introduction 

We have .seen in the last chapter that there remains a few anomalous features such as the 
large linewidth, significant asymmetry in the Raman lineshape and substantial intensity ob- 
served in the classically forbidden Ai g geometry could not be explained within the existing 
theoretical framework. (1 6] Several new ideas have been put recently to resolve these issues 
which include the exchange disorder caused by zero point lattice fluctuations, [7] and resonant 
Hainan scattering in which the incident photon energy is comparable to the energy gap in 
cuprates[S]. 

When lattice distortions, caused by quantum and thermal fluctuations, are taken into 
account, in the adiabatic or Born-Oppenheimer approximation, the hopping terms t t] , and 
hence the exchange couplings Jp, which depend upon the instantaneous position of ions, 
must therefore include random terms. This approximation should be valid for high-energy 
magnou excitations in cuprates where the rnagnon energies of order 2 J ~ 2000cm -1 are much 
greater than the Debye energy hui® rsJ 340cm -1 . The exchange disorder caused by zero-point 
lattice vibration was recently taken into account within a nearest-neighbor (NN) Heisen- 
berg model where the exchange energy J + SJij at each bond includes random terms. [7] For a 
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Gaussian distribution of the random terms, this model was studied using the Quantum Monte 
Carlo (QM< ) technique and Raman intensities were obtained numerically in different scat- 
tering symmetries. Satisfactory fitting was obtained with the experimental Raman intensity 
lineshape for a Gaussian variance of a ~ 0.4. However, no clear insight emerges from this 
QMC calculation regarding the explicit effects of exchange disorder on (i) low-energy, long- 
wavelength magnon modes, (ii) the magnon velocity, and (iii) the high-energy magnon modes 
with energy 2,/ which are mainly responsible for the Raman scattering. Ascertaining the 
exchange-disorder-induced renormalization of the magnon energy scale is of vital importance, 
however, as it is from comparison of this scale with experiments (such as neutron scattering, 
Raman scattering, magnetization versus temperature etc.) that a reliable value of J for the 
cuprates is extracted. 

In this chapter we therefore examine these questions in the context of an antiferromagnet 
(AF) with off-diagonal disorder. We consider the following Hubbard Hamiltonian on a square 
lattice’, with randomness in the hopping terms, and with a filling of one fermion per site, so 
that an antiferromngnetie insulating state is obtained. Generalization to three dimension and 
to other bipartite lattices is straightforward. Random terms Stij are included in NN hopping 
terms, and are chosen independently for each NN pair of sites from a Gaussian distribution, 


K ■- — ^ {t + 5tij){a\ (r aj a + h.c.) + U fiit^ii, (3-1) 

(«»."■ ‘ 


PiStJt) = 


1 

\/2na 


exp 


2 a 


(3.2) 


The distribution width \Ja measures the disorder in hopping. Since for weak disorder 
~ 2 therefore the distribution widths y/a] and y/F t for random exchange-energy 
terms and hopping terms should be related via crj — 4 oq. An estimate for crj has been 
obtained(7] using crj « {()lnJ/dr){{5r)/r), where the relative amplitude of the zero-point 
motion (hr)/r was estimated as about 5%. The dependence of the exchange energy J on the 
in-plane ( ud ’u distance r is taken to be of the form J ~ r"°. The exponent a for cuprates 
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has been obtained using high-pressure studies[9] in the range a = 4 ± 2, whereas a stronger 
dependence (a as 6.8 ± 0.8) has been reported by others.[10] In their QMC calculations Nori 


ft at. have considered the cases oj = 0.3, 0.4 and 0.5. 

We have examined — both perturbatively and numerically — effects of hopping disorder 
on magnon energies, their wavefunctions, and on the density of states (DOS). In both schemes 
the magnon modes are obtained at the level of Random Phase Approximation (RPA). At this 
h'vel the magnon propagator has the form (x“ + (w)] = [x°( w )]/l “ ^[x°( w )]> where [x°( w )] i 8 
the rerot h -order antiparallel-spin particle-hole propagator evaluated in the broken-symmetry 
AF state. Eigensolutions of the [xV)l matrix yield the magnon energies and wavefunctions. 
fr. the perturbative analysis we obtain the disorder-induced perturbation to [x 0 ( w )]» 
obtain resulting corrections to its eigenvalues, which then yield the renormalized magnon 
energies. For analytical convenience we have considered the strong-correlation limit as the 
perturbative analysis is particularly simple in this limit. We have also numerically obtained 

\ *)I using the self-consistent eigenstates of the HF Hamiltonian in the AF state with the 

r ,md om hopping terms included from the beginning. This method can be applied for arbitrary 


interaction strength U/t. 

\W find that there is strong renormalization in energy of high-energy magnon modes due 
, u - t cooperative effect arising from local correlations in hopping disorder, which results in 
appreciable magnon DOS well above the maximum energy 2 J for the pure system, even for 
r small as 0.1. This result is significant because it not only clearly shows how the two- 
rr.»«uon Raman scattering intensity, which involves the one-magnon DOS, can extend well 
hrvo nd the energy 4 J when hopping disorder is present, but it also yields an insight into why 
R.rix.-au scattering is so sensitive to the tiny zero-point fluctuations. We also show that these 
•; ir a* - energy magnon modes at the upper end of the spectrum are strongly localized in regions 
lattice where the locally-averaged hopping strength is maximum, i.e. at sites where 
^ if./s connecting to the site have maximum average. We have also obtained the energy 
Z. > realisation of the low-energy, long-wavelength modes, both to first order and second 
-r in BUilt. We find that the first-order correction to magnon energies are of the random- 
t ype, which decrease like IfVN with increasing system size, but can be significant 



for finite-size lattices Tk» i i 

ne second-order correction, on configuration averaging, yields a 

ntum independent multiplicative renormalization of the magnon energy by a factor 1 +a, 
- ( + cr)yJ\~Z 7 ^ Thus the Goldstone mode is preserved, and the low-energy, 
g avelength magnons are weakly perturbed by off-diagonal, hopping disorder. This is 
contrast to the case of potential scattering (diagonal disorder) where in two dimensions 
g Iar corrections to low-energy magnon modes were obtained[ll]. 

3.2 Perturbative Analysis 

obtain the magnon energies and wavefunctions, we evaluate the transverse spin fluctu- 
ation (magnon) propagator *- + K r't') = <AF|S-(ri)S+(r'*')l AF) which, at the RPA level, 
is given by (\ (w)j/l - £/(x°(u;)] after Fourier transformation to frequency space. Here 
[\ (^ )] is the antiparallel-spin particle-hole propagator evaluated in the broken-symmetry, 
self consistent state, with matrix elements given by, [x°(^)]ii — i /(^ 72 ^)G r J J (u; , )<?j - ,(aj , -a?). 
The magnon energies are then given by the poles in the magnon propagator, 1 — UX(ui) = 0 , 
where A(w) is the eigenvalue of the [x°(w)] matrix. 

In the perturbative technique, the disorder-induced perturbation [£x°] to the [x°(u>)J mat- 
rix is obtained diagrammatically in powers of SUj/t. Resulting corrections to the eigenvalues 
^ ^ en yield the renormalization in magnon energies. Diagrams which contribute 
to (V) to first order in Sty are shown in Fig. 3.1 (a). The upper and lower lines are re- 
spectively the t-spin and j-spin propagators with frequencies u/ and u/ - u. In the diagonal 
matrix element [fy 0 ];., k is summed over all NN sites of i. Other diagrams are obtained by 
exchanging the * and k vertices of the dashed line and by putting the dashed line in the lower 
(j-spin) propagator etc. In the strong-correlation limit we obtain, 


= ~2A3 — 2A 3 ^* j ^’ (3-3) 

where 2A = mU is the Hubbard gap characterizing the AF insulating state, and only terms 
upto order (t 2 /A 3 ) have been retained, appropriate to the strong-correlation limit. We notice 
that the sum of all the NN matrix elements [Sx%j I s precisely the diagonal matrix element 
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[<J\ ]i(. An immediate consequence of this correlation is that the Goldstone mode is preserved 
and that generally the effective scattering of low-energy, long-wavelength modes is weak. To 
begin with we consider the first-order correction to the eigenvalue A(w) due to the perturbation 
[£\ 0 ]. For this we require the eigensolutions of [x°(uj)] in the pure AF state, which can be 
labeled by momentum due to translational symmetry within the two-sublattice basis. We 
first consider the case of low-energy, long-wavelength modes for which 1 and w/2 J 1. 
I'pto order q 2 , (u/'U) 2 and in the strong-correlation limit, the eigenvector and eigenvalue of 
(\°(u;)] for the pure AF are given by, [12] 



\ A + u i 2J 'i 

) 


(3.4) 


>(o> = I _ Jl _ I (!£.)* 

5 U A 3 [4 2 \2JJ 



(3.5) 


• • * • ■ * 
^ 1 | ^ 1 



Figure 3.1: Diagrams contributing at the (a) first-order level in St (dashed lines) and (b) 
second-order level to the diagonal matrix element [<$x°]» nearest-neighbor matrix 

element [d* 0 ]# 

Therefore for the first-order correction we obtain, 
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( r * + r ») (jj) + + r v ql 


(3.6) 


where r x = £ £• 8t {j /t (fj = r\- + ax) 

and r y = Stij/t (fj = r { + ay) (3.7) 

are the lattice averages of the random hopping terms Stij taken in the x and y directions 
respectively. Significantly the order u terms cancel exactly for arbitrary distribution of 8t,j, 
and their cancellation follows from the fact that the A-and B-sublattice sums of the diagonal 
terms are identical, i.e. £, € a[^X°]« = E.'eB^X 0 ].*- 

Now, r r and r y will vanish for an infinite system if the random hopping terms Stij are 
distributed symmetrically about zero, and therefore in the infinite-size limit, there will be no 
first-order correction to long-wavelength magnon energies. However, for a finite-size system 
the averages r x and r y are like random-walk averages and will scale like l/y/N. This is 
because the terms £, SUj/t in r x and r y are sums of N Gaussian random variables, and 
therefore will themselves have a Gaussian distribution of width \fNa. 

Solving for w from the equation, 1 - U{Xf ] + 5AJ 1 )) = 0 yields the renormalized magnon 
energies for long-wavelength modes: 


w, = y/2 J[ql{ 1 + 3r* + r y ) + gj(l + 3r y + r*) 1/2 ] (3.8) 

We now examine the effects of hopping disorder to first order on the high-energy magnon 
modes with energy « 2 J. In the pure AF the maximum magnon energy is 2 J and these 
magnon modes have non-zero amplitudes only on sites of one sublattice, corresponding to 
creating spin deviation on these sites. The energy cost due to a spin deviation of this kind is 
therefore 4(J/2) = 2 J, where J/2 is the bond strength and 4 is the number of broken bonds. 
Now, when hopping disorder is present the highest-energy magnon mode will correspond to 
creating a spin deviation on that site where the sum of the NN bonds is maximum. This 


will clearly occur in a region of the lattice where the hopping disorder terms SUj neighboring 
a site have the maximum average. In fact, generally the high-energy magnon modes with 
energy > 2 J will be localized around such sites across the lattice. 

Interestingly, in such regions where the locally averaged hopping is significantly higher 
than the bulk average f, the diminished U/t ratio leads to a lowering of the local staggered 
magnetization m(r), and hence of the local charge gap 2A(r) = m(r)U. Fig. (3.2) shows the 
magnetization m(r) calculated for a — 0.1 and U/t = 10 for a 10 x 10 system. As mentioned 
above we can see a local lowering of the charge gap at certain regions where the magnon 
wavefunction is found to have the maximum amplitude. Thus the localization of the high- 
energy magnon modes in such regions is suggestive of the high-energy magnons acquiring 
a charge due to hopping disorder and getting trapped in the local depressions of the charge 
gap. This can be seen formally by examining the self-energy correction, Ey = U 2 5x*ij, to the 
magnon modes due to hopping disorder, particularly the diagonal terms which will act like 
potential terms for the magnons. From Eqn. (3.3) we obtain, 

= 

Now *£.& &U,i+s/t the sum of z (coordination number) Gaussian variables and will therefore 
itself have a Gaussian distribution with an effective second moment of Zcr. Thus the local self- 
energy S„ has a scale of J and is randomly distributed according to a Gaussian distribution 
with a much bigger width \ZZcr. The high-energy magnons thus effectively experience a 
random potential and therefore get trapped in the local potential depressions. In fact, E« 
precisely yields the extra magnon energy for the localized high-energy modes, and therefore 
we can conclude that these modes will extend well beyond 2 J in energy, and that the excess in 
energy will have a Gaussian distribution of width VZa- This result has a significant impact 
on the one-inagnon density of states, indicating that there will be a Gaussian tail m the DOS 
extending well beyond energy 2 J. Indeed, this is evident from Fig. (3.3) showing the magnon 
DOS obtained numerically which is discussed in the following section. 

We have thus shown that while the low-energy, long-wavelength magnon modes are very 


s 


0ti,i+5 


(3.9) 
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Figure 3.2: The amplitude for the highest magnon-mode (below) and the magnetization profile 
(top) obtained from numerically is shown for a 10 x 10 lattice for U/t = 10 in presence of 

disorder ( rr — 0.1). 

weakly affected by hopping disorder, on the other hand the high-energy magnon modes are 
strongly alieeted because of a cooperative effect due to local correlations in hopping disorder. 
There is therefore a strong asymmetry in the energy shifts of low-energy and high-energy 
magnon modes, which is also evident from the magnon DOS shown in Fig. (3.3). Now the 
asymmetric iineshape in the two-magnon Raman scattering in cuprate antiferromagnets has 
been a puzzling feature. If indeed it is the randomness in hopping due to zero-point motion 
which is responsible for this feature, then we have shown that this asymmetry is basically 
arising from the asymmetric energy renormalization of high-energy and low-energy magnons 
due to hopping disorder. Also, the strong energy renormalization of high-energy magnons 
due to this cooperative effect of hopping disorder leads to an insight into the extraordinary 
sensitivity of the Raman scattering process on zero-point motion of the lattice. This is the 
central result of this chapter. 

We have seen that the first-order correction to [x°M] yields vanishing contribution to the 
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iMj’irn- 3.3: Tlu* configuration-averaged magnon density of states obtained numerically for 
a 10 >: 10 hopping disordered system with a = 0.1 and 0.3, in the strong-correlation limit 

(1'i't 101. 


low energy, long- wavelength magnon energies in the limit of infinite system size. Therefore it 
is necessary to consider second-order correction to [x°(ce)], for which the relevant diagrams 
are shown in Fig. 3.1(h). Humming all the processes, we obtain after configuration averaging, 



hVw)|=il-(l + cr)|l l_2l] (3-11) 

where 2.1' s 2./(l + a) is the renormalized magnon-energy scale. Therefore, all magnon 
energies are simply renormalized by a momentum-independent multiplicative factor (1 + a), 

and we have, 


w 9 = 2J(1 + <7)(1— t, 2 ) 1/2 


(3.12) 
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3.3 Numerical analysis 


We now briefly describe a numerical method for obtaining the magnon energies and wave 
functions. Details of this method have been discussed earlier for the impurity-doped AF 
with on-site potential disorder. [13] The idea is to self-consistently obtain the AF state for . 
the hopping-disordered system within the Hartree-Fock approximation, and then determine 
the magnon eigensolutions from the transverse spin propagator, evaluated in the broken- 
svmmet rv state. Since wo are concerned here with an AF insulating state for the one-fermion- 
per-site system, the II ¥■ plus-fluctuations approach is expected to be quantitatively accurate 
even in the strong correlation limit, as has been shown to be the case for the pure system. [14] 
Within the HF approximation we construct the Hamiltonian matrix on a two-dimensional 
scpiat-e lat.t ice with periodic boundary conditions. The random hopping terms 6tij are obtained 
for each pair of NN sites using a Gaussian random number generator. In the site basis 
matrix elements of the Hamiltonian for spin a are given by (i\H„\i) = U{n iW ) and {i\H„\j) = 
..{t {fit,,) if" i and j are nearest neighbors. Starting with some initial configuration for the spin 
densities, the Hamiltonian matrix is constructed and diagonalized to yield the eigensolutions 
| /G.o,,,}. From these the spin densities are re-evaluated, the H-matrix is updated, and the 
ptoeedure iterated until self-consistency is achieved. These eigenstates obtained numerically 
of* the self consistent AF state are used to construct the [x°Hl matrix as described earlier.[13] 
It is convenient to express the transverse spin fluctuation propagator at RPA level in 
terms of the eigensolutions {<i*\{w), A(to)}, of the [xo( w )] matrix. 


(x~»] = E 




(3.13) 


A x - U A(w) 

In this form we have a propagator representation for the magnon, with the magnon wave- 
f, given l,y Ur eigenvector |*(r)> and magnon energies obtained from the pole, 
,, o. 'll,,, energy w„ of the n-th magnon mode is obtained by solving, 1— CA n (to„) = 
(I for the appropriate (n-th front the top) eigenvalue A.M- The root of the equation, 
A„U,! =. i/f is determined by obtaining A„(w) for closely spaced values of u, on both 
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sides of the root, and finally linearly extrapolating between them. Suppose (A‘,ud) and 
( A n'U-fJ are two sets of values for two energies very close to, and on either side of the root. 


then the root u.’ n is determined from, 

A(u.’„) = - 


>1 ,%-K 


(3.14) 


Tin' whole procedure therefore involves obtaining the self-consistent state for a given realiza- 
tion of the hopping-disordered system, followed by constructing and diagonalizing the [ x °)( u .-j] 
matrix, and finally from the set of eigenvalues A n (w) the magnon energies are obtained by 
intt i point ion as dest ida d above, This scheme for obtaining the magnon energ i es is then 

repeated for several configurations. 

We find that the highest magnon energies always extend well beyond the maximum 
energy 2.7 for the pnn' case. Even for a as low as 0.1, the highest magnon energy obtained is 
1.3d in units of t whereas 2.1 ~ 0.8 in units of t (for U/t = 10). This is in agreement with the 
perturbative analysis finding that the high-energy magnon modes are strongly renormalized 
due to the correlated effect of hopping disorder. 

In Table 3.1 we show the variation with a of this hopping-disorder-induced increase in 
the magnon energy for the highest magnon magnon mode, and compare with the perturb- 
ative result derived earlier (Eqn. (3.9)]. For a given a the highest magnon mode energy 7 
is configuration averaged over 25 configurations of the hopping-disordered system to yield 
l\n. The difference with the maximum magnon energy for the pure case yields the excess 
due to hopping disorder. This difference is compared with the perturbative result for the 
magnon self energy given in Eqn. (3.9) for high-energy magnon modes with energy near 2 J. 
For a given realization of the random hopping terms, the nearest neighbour sum Yls$tij/t is 
performed for ail lattice sites, yielding the local self-energies an. From this the maximum is 
picket! out, and this is averaged over 25 configurations, yielding the configuration-averaged 
maximum self energy correction, a^. That the excess is magnon energy for the highest mode 
matches well with the maximum of the local self-energy correction strongly supports that the 
niagnott modes involves spin-deviation at those sites where the locally averaged 
hopping terms are the maximum. This implies that the upper end of the magnon-energy 
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spectrum is entirely determined by local correlations in random hopping terms. In as much 
as this part of the magnon spectrum is most important in the two-magnon scattering, fea- 
tures such as the strong asymmetry can thus be traced to these local correlations in hopping 

disorder. 


'fable 3.1 : Yariaf ion of the configuration-averaged maximum magnon energy tZ^f with hopping 
disorder strength a. Also shown are the configuration-averaged increase in magnon energy 
yy n ... for the highest magnon mode relative to the pure-case energy te^, and compar- 
ison with the configuration-averaged maximum self-energy correction S m evaluated from Eq. 

(3.9b 


(T 

u " J m / J 

‘J/J 

E m /J 

0.0 

1.89 

0 

0 

0.1 

3.05 

1.16 

1.12 

0.2 

3.44 

1.55 

1.58 

0.3 

3.95 

2.06 

1.93 


A convenient way to exhibit the strong renormalization of high-energy magnons is via 

the magnon density of stall's: 


1 1 

•'V) = I'SzTirZJ 


(3.15) 


T)~+0 7 T Nn n (id ~ U n ) 2 + 7] 2 

where _o..'s denote the magnon energies and r; is chosen to be of the order of average level 
spacing. ,V„ represents the total number of modes obtained from the 50 different configur- 
ations considered. The results for the magnon DOS are plotted in Fig. (3.3) for different 
values of a. Also shown for comparison is the DOS for the pure system which extends upto 
energy - 2,7. It is observed that the DOS for the hopping-disordered system extends well 
beyond the ma Amum energy 2.7, though the position of the peak remains roughly the same. 

We have also examined the magnon wavefunction for modes at the high-energy end of 
the son t mm. and find that indeed these modes are strongly localized in certain regions 
„f til,- lattice. « expected from the analysis discussed earlier. The magnon amplitudes 
are contained in the eigenvectors |*(«,)> * “>« Wl matrix « wta ** ed * 



energies. A useful way to exhibit the magnon modes on the lattice is via arrows with size 
proportional to the local magnon amplitude <£) v , and orientation proportional to the rotation 
angle 0, = sin“ l (<?\/5l). where Si is the local magnetization. [15] The zero-energy (Goldstone) 
mode in this representation is characterized by identical rotation angles, despite the local 
magnetization and magnon amplitudes having different magnitudes on the lattice sites due to 
disorder. We have also confirmed that this localization of high-energy magnon modes occurs 
in those regions of the lattice where the locally-averaged hopping is significantly higher than 
the bulk average, so that the locally-averaged Ujt ratio and the local magnetization are low. 

In conclusion, we have studied the antiferromagnet with hopping disorder, both per- 
turbatively ami also using a numerical method, in the strong correlation limit. We have 
shown that while the low-energy, long- wavelength magnon modes are only weakly affected, 
the high-energy modes are strongly affected because of a cooperative effect arising from local 
correlations in hopping disorder. There is therefore a strong asymmetry in the energy shifts of 
the low-energy and high-energy magnon modes. If indeed zero-point motion leads to hopping 
disorder in cuprates, then this work provides a qualitative understanding of the anomalous 
linewidth ami asymmetric lineshape seen in the two-magnon Raman scattering in cuprate 
antiferromagnet*. Also, the strong energy renormalization of high-energy magnons due to the 
cooperative effect of hopping disorder leads to an insight into the extraordinary sensitivity 
of the Raman scattering process on the zero-point motion of the lattice. 



Appendix A 

In this Appendix wo discuss the evaluation of the matrix element in the strong correl- 

ation limit wherein terms only upto order have been retained. From Fig. 3.1(a) one can 
write. 


M* = i / (3.16) 

where G' s represent the host Green’s functions, 0,© refer to the quasiparticle bands and the 
sum over k represents the sum over all z nearest neighbours. Substituting the expressions 
for the Green's functions appropriate upto order tj A one obtains, 


(£\o(**)]„ ® 1 ] 2 ^ 


-t/2A tj 2A 1 ( 1 \ 

u) 1 — A® — iff w' — A® + irj J —ui — A® + irjJ 


(3.17) 


where A' '/S r ~ - A/ + A. 


The frequency integral yields, 

W4 = -S3Et < U8 > 

Identical result is obtained corresponding to the dashed line (Silt) being put in the lower (!- 
spin) propgator. Hence summing up the contribution from both these processes we obtain, 


r ^ 

“ “2 A 3 ■ 


SUk 

t 


Similarly the off-diagonal contribution ($Xo] ti can be shown to be equal to, 


[$XoM]y 


£ % 

2A 3 t 


(3.19) 


( 3 . 20 ) 
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Chapter 4 


Hopping-disorder-induced effects upon 
Raman Scattering spectrum in an 
Antiferromagnet 

4.1 Introduction 

la Chapter 3 we have discussed the hopping-disorder-induced effects upon the magnon spec- 
trum in an antiferromagnet. It has been shown explicitly that while the low-energy magnon 
modes are negligibly affected by disorder, there is a substantial effect upon the high energy 
modes. It was also mentioned there that such asymmetric magnon energy renormalization 
can provide a clue to the considerable asymmetry in the two-magnon Raman lineshape. Here 
we are going to investigate the effects of hopping disorder upon the two-magnon Raman 
scattering intensity in an antiferromagnet. 

We have seen earlier that in Nori et aV s[l] approach the hopping disorder is looked upon 
;is an interaction between the lattice excitation (i.e. phonons) and the spin excitations (i.e. 
magnons) which induces changes SJ ih in the exchange integral J. An estimate of the mean 
square deviation (C-j 1 ) 2 ) is made and to make a comparison with the experimental data the 
<:or„. S po,.,lh.g value is considered to be as high as 0.5. It has been suspected[2] that such a 
large variation in J appears physically unreasonable and questioned that whether more mod- 
erate values of (^) are sufficient to achieve agreement with experiments. Essentially what 
one obtains is that the coupling of phonons to the spin excitations leads to the broadening of 


*7£ 
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a two-rnagnon peak, the asymmetry of the lineshape about its maximum towards the higher 
eiu'ig} regime and appearance of scattering intensity in forbidden geometries like A lg , B 2 g etc. 

hai lim , we ha\e investigated, [3] the effect of lattice fluctuation induced hopping disorder 
in the context of Mott-IIubbard AF where the hopping term tij include random terms. The 
Hubbard Hamiltonian in this case can be written as, 


// = 


- S (* + + h.c.) + {7]T 

i 


(4.1) 


{Tj = 7?Ta exp 


-{Halt? 

2 a 


(4.2) 


where the' random terms Sti 3 ' s for each bond are chosen independently from a gaussian dis- 
tribution and the distribution width yfa measures the strength of disorder. There we have 
examined both perturbatively and using a numerical method - effects of hopping disorder on 
mugnon energies, their wavefunctions and on the density of states (DOS). In both schemes, 
the spin- wave energies are obtained at the Random Phase Approximation (RPA) level. We 
have obtained a strong renormalization in energy of the high energy magnon modes, arising 
from locally correlated hopping, which results in appreciable one-magnon DOS well beyond 
the maximum spin-wave energy 2 J for the pure system. This result has got important im- 
plications in explaining the experimentally observed Raman spectrum. 


To obtain the Raman scattering intensity we need to consider the imaginary part of the 
full two magnon propagator which in turn requires the knowledge of spin-wave energies. For 
a finite 1 size system these energies are obtained numerically at the RPA level in presence of 
disorder using a numerical method discussed in detail in Chapter 3. Following a prescription 
of the mode assignments discussed in section 4.3, it is possible to obtain the symmetry factor 
required for the B lg scattering geometry. It is observed that even for a sufficiently small 
strength of disorder, a = 0.01, (much smaller than the values considered by Nori et a/.[lj) 
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though the Raman spectrum on the lower energy side is weakly affected, there is a substantial 
( hange in the highei fietjiu nc\ icgimc and a good agreement is obtained with the experimental 
data. Ibis confirms om eailiet statement [3] that if indeed the randomness in hopping due 
to zero- point lattice fluctuation is playing a key-role in explaining the long standing puzzling 
feature of tin* Raman lineshape regarding the asymmetry, then it essentially arises from the 
highly asymmetric magnon-energy renormalization. 

4.2 Configuration averaging of the two-magnon propag- 
ator 

The Kleury-I.oudon (FL) Hamiltonian which represents the interaction between the photon 
and the spin pairs can he written as, 

U It = A £(E in <J)-(E sc i)S(r).S(r + 5) (4.3) 

r ,5 

Where the E;,ic and E,- are the incident and the scattered electric field vectors, 6 is the unit 
vector connecting nearest neighbour (NN) sites of opposite sublattices and A is a constant. 
The sum over r ensures that the spin excitations have net zero momentum as the incoming 
photon has essentially zero momentum. For calculation of the Raman intensity one needs to 

consider the correlation function of the spin pair operator, 

P$ = S(r).S(r + 5) (4.4) 


For the Bi„ symmetry we have, 



4.2 Configuration averaging of the two-magnon propagator 


79 


Emc ^ x y , 
E iC ~ x + y. 

So that the Hamiltonian assumes the form, 


Hr ~ S(r).{S(r + px) — S(r 4- py)} 

r,p=±l 


(4.5) 


S(r) 


+ 


S( ri ) S( ri ) 


ri 


S(t’) 


r + 5 


•- 

+ . 


S„'(r) 

B Ig 


ri+jLL 


S B (rj) 

Big 

S (h) 

Aig 


S B + , (ri) 

S (r,) 

Aig 


r’ + 5’ 


s' <y) 

Big 


1.1: Interaction between two magnons is shown where the bold lines represent the 
two magnons propagating from neighbouring sites r and r + 5 to r 1 and r' + d'.The wavy 
lines denotes the interaction between the magnons at sites r x and r x + /r where fx connects 
t„ the nearest neighbouring sites. The creation and destruction of magnons at these points 
are shown by 5'+ and S~ where the subscripts Big and A Xg correspond to different scattering 

geometries defined in the text. 


Saw lot us consider the diagram representing (Fig. (4.1)) the interaction between the 
magnons in real space. The interaction between the magnons excited at sites r x and n+/x 

can be written as, 
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V,„t - -J Y, $ { r i)S + (r 1 )S (r x + iu)S + (r 1 + y.) (4.6) 

1*1 

where fi connects to NX sites of the opposite sublattice. The form of the interaction indicates 
that the propagators on either sides of the interaction line are coupled as y can connect both 
in the i‘ and ?/ direction. But. they can be decoupled if one rewrites the interaction term as 

(•see Fig. (4.1)), 


U„--./E 5 '(r 1 )S + (r l )S 5 1 (r,)SJ,(r 1 ) + S^rOS^r,) 

+5o 3 (r 1 )S^ 3 (r 1 ) + So 4 (ri)5<? 4 (r!) 
where the 4 operators are defined in the following way, 

s oM) = Iy,{ s± ( t + p*) - s± ( r + py)} 


(4.7) 


$0,(r) = ^ + s± ( r + py)} 

~ p 

S'oJ r) = \ Yl{P S± ( r + - P S± ( r + PF)} 

w p 

S ( \( r) = ^£{P 5± (r + p£) + pS ± (r + py)} 

- p 

where the 4-sign and the -sign in the superscript refers to raising and lowering operators 
respectively. It can be seen that these operators correspond to the the different irreducible 
representations predated in Appendix C of Chapter 2. The first two operators correspond 
to B 1r and A,* symmetry factors respectively. Hence explicitly we write them as S Blg and 
St for our convenience. In this new notation the interaction term looks like, 
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U» = -./E.f-(r l )if + (r 1 )[5i li (r 1 )Sj, i (r 1 ) + SXJr.JS+Jr,) (4.8) 

+S5,(ri)S$,(r l ) + S5 4 (r,)SJ ( (r t ) 

Having written the interaction term we proceed to define the spin-wave propagators for 
different geometries. In Fig. (4.1) the spin raising and lowering operators at r and r -f 5 
are denoted hy >~(r) and N Big (r) and those at at r' and r' + S' are denoted by S~( r') and 
■Sii- ( r/ )- f ids immediately suggests that 5 and S' which connect NN pairs are fixed in space 
and we can write the spin-wave propagators as, 

X - + (r,r') = (V>d5-(r)5 + (r / )|0 G ) (4.9) 

corresponding to the upper line and. 


\B~( r ’ r ') = {V’G|S^ g (r)5 Blg (r , )|^G) (4-10) 

corresponding to the lower one in Fig. (4.1). Now the time-ordered two-magnon propagator 

in the B.„_ geometry can be defined as, 


( r ,«)S + (V,t')SS„(r',«')Wo) (4.11) 


I.'or caliuiat ing the Raman intensity first consider the non-interacting limit of the two-magnon 

propagator. 



I'he configuration averaged two-magnon Raman propagator is given by, 


(4.12) 



(4.13) 
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Wo can decouple the two propagators if we neglect the corrections arising due to exchange 
between the mngnons. This allows us to write, 


G\ 


+ ( I b r '> w i) XB lg (r,r',u; -wi) 
r,r' 


(4.14) 


Having done the configuration averaging of G°(u) one can proceed to calculate the two- 
magnon Hainan intensity which is obtained from imaginary part of G(u) where, 


G o(cu) 

1 + J Gq(u) 


(4.15) 


We can neglect the vertex corrections since they appear only at the second order level and 
hence proportional to the fourth power of the disorder strength he. (y) 4 (in reference [4] 
this result is shown explicitly for the on-site (diagonal) disorder case and calculation of the 
vertex correction is discussed elaborately in presence of on-site disorder in Chapter 5). So 
for weak disorder, the correction will be insignificant. This enables to write, 


G_ o(^) 

1 + JGju) 


(4.16) 


We have shown it, in Appendix A that the B lg symmetry factor viz. (cosq x - cosq y ) 2 
, pondine to the pure case is obtained by Fourier transforming the non-interacting two- 
magnon propagator ('4, . This symmetry factor is important for calculating the two-magnon 
Hainan pe< trum. It can also be noted that due to the form of the interaction i.e. presence of 
both B u . and A,„ terms in V mt , off-diagonal combinations of different scattering geometries in 
f are possible. In that cast* the symmetry factor will contain contributions from different 
geometries. We can see it clearly if one considers the term, (S'^ g (r)5X lf (r , )) ) which when 
Fourier transformed yields, 


^(cos<? x — cos <? y )(cos q x + cos q y )x + (4-17) 

>1 

But this symmetry term can be shown to be identically equal to zero for all frequencies be- 
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cause of the symmetry present in the magnetic Brillouin Zone. Hence the only terms which 
survive are the ones that are diagonal in scattering geometries [5]. 

4.3 Mode Assignments 

Hen' we want to establish a one-to-one correspondence between the pure system and the one 
in presence of hopping disorder. This we believe to be valid if the strength of disorder is 
weak. For doing this we calculate the magnon modes numerically for a 10 x 10 system in 
absence oi disorder in the strong-coupling limit (UJt = 200) . This procedure of obtaining 
the magnon modes numerically has been discussed in Chapter 3. These modes thus obtained 
are assigned a pair of moth' numbers n x and n y which are related to the wavevectors q x and 
<!,, via the relation q x / y = qrn x / y in the translationally invariant (pure) case. The numbers 
;>j. and n y are chosen such that the corresponding q x and q y are restricted to the upper-half 
(first two quadrants) of the Brillouin Zone (BZ) where q x is allowed to have both positive 
and negative values (ranging from — tt to +7r) whereas q y is assigned only positive values 
(from 0 to rr ) . Tin* other half of the BZ does not need to be considered as the q x and 
q, t values (or equivalently the n r and n y ) lying within the region specified (see Fig. (4.2)) 
are sufiieieut to label ail the modes. The Goldstone mode or the lowest energy mode is 
assigned with the mode numbers n x = 0, n y = 0, the second collective mode with n x = 1, 
n, t ■ » (or equivalently n r = 0, n y = 1), the third one with n x = 1, n y = 1 and so on. 
Having done with the assigning scheme it is now possible to form the symmetry factor for the 
B t , scattering geometry which is = [cos(^)n x - cos(^)n y ] 2 . As mentioned earlier this 
svmmet rv factor is crucial for calculating the two-magnon Raman intensity. The consistency 
of the mode assigning scheme is checked by re-evaluating the magnon energies using the 
mode numbers via the relation a i n = 2j/T^7 2 where 7n = |[cos(f )n x + cos(f )n„] and 
./ r. s p 0.0 2 (i.e. Vjt = 200). The energies thus obtained show excellent agreement with 
the ones calculated from the numerical method. The same mode assignments are expected 
to be carried over for the intermediate coupling regime. We show the comparison between 



4.3 Mode Assignments 


84 


the energies calculated numerically and the ones evaluated by using the analytical expression 
relevant to the intermediate- £7 regime (U/t = 10) in Table 4.1 where we have, 

1/2 

(4.18) 


2^ U-7,0 


-{6 + 3cos 


n x cos 



1 .1 2: The region in the Brillouin Zone (BZ) is shown where the dots denote the mode 

and n, chosen so as to assign labels to to the mode energies calculated. A few of . 
mud,, numbers are indicated. The first two quadrants are required for assigning numbers 

to all the modes. 


We establish the correspondence 


between the pure and the disordered system by assuming 
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the same assigning scheme continues to work for the disordered system. We show the mode 
energies for the disordered system (a = 0.1) for a particular configuration in Table 4.2. We 
also notice that the degeneracy of the modes is lifted due to the presence of disorder and the 
low eneig\ modes aie not significantly affected whereas there is a considerable effect towards 

the higher energy side. 

lable 4.1: I he mode energies are obtained using the expression relevant to the intermediate- 
l regime (l /t ~ 10 or A // = 5) for a 10 x 10 system where n x / y = (2k )q x / y are the 
mode numbers specified in the leftmost column. Also we have shown the energies calculated 
numerically using the exact eigenstate method for comparison. The symmetry factors for 
both Bi ;> and .ii 3 scattering geometries are obtained using the mode numbers. 


Mi )« U % Nos. ( iij , n 1; ) uj n 

(numerical) 

uj n (formula) 

K(B 1.) 

4&Ai a ) 

C-o.O) 

0.7548 

0.7505 

4.000 

0.000 

(4.01 (0.4) (-4.0) (-5,1) 

0.7438 

0.7408 

2.618 

0.000 

(4.1) (1.4) (4,1) (-1.4) 

0.7396 

0.7356 

0.382 

0.000 

',o.2i 1 2,5) ( o,2) ( 2,3) 

0.7150 

0.6931 

3.272 

0.037 

(3.U (1.3) ("3,1 ) (-1,3) 

14,21 (2,4) ( -4,2) (-2,4) 

0.7030 

0.6927 

1.250 

0.250 

1 3.0 1 i 0.3 i I 3.0) (-5.2) 

0.6882 

0.6793 

0.000 

0.382 

'( 2 . 2 ) i:i,:V)'('-'2.2)'"('-3".3'r 

0.6765 

0.6696 

1.713 

0.478 

(2,iV(i,2f("’2:i) vn— 

(4.3) (3,4) (4.3) (-3,4) 

0.5945 

0.5826 

0.250 

1.250 

1 2,0 1 1 0,2) ( 2,0) (5,3) 

0.5469 

0.5458 

0.477 

1.714 

( i.h (4.4) (1.1) {-4,4) 

0.4199 

0.4102 

0.000 

2.618 

(1,0) (0.1) (1,0) (-5,4) 

0.3061 

0.2987 

0.036 

3.273 

(«K0) 

0.0000 

0.0000 

0.000 

4.000 
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lable 4.2: Here the assignment of the mode numbers n x and n y is shown for the mode energies 
calculated numerically for both pure and disordered systems (cr = 0.1) of system size 10 x 10 

and for l'/t = 10. 


Mode NoS. (lljr.Uy) 

ui n (pure case) 

ui n (disordered case) 

(-5,0) 

0.7548 

1.3006 

(4,0) (0.4) (-4.0) (-5,1) 

0.7438 

1.2222, 1.2024, 1.1960, 1.1346 

(4.1) ton (-4,1) (-1,4) 

0.7396 

1.0578, 1.0178, 1.0137, 0.9955 

0.0 0-;n 1 4.2) 1-2,4) 

0.7150 

0.8557, 0.8414, 0.8291, 0.7934 

(4.1) 11.4) ( 4,1) (-1,4) 

0.7030 

0.7854, 0.7762, 0.7698, 0.7641 

(4.2) (2,4) (-4,2) (-2.4) 


0.7228, 0.7018, 0.6935, 0.6624 

"(4.0 ! (0.4) 1-4.0) (-5,2) 

0.6882 

0.6614, 0.6554, 0.6364, 0.6197 

72,2)1 :U)1"2.2) (-4,4) 

0.6765 

0.6020, 0.5773, 0.5668, 0.5411 

lllTTcii) (-2,1) 1-1.2) 

0.5945 

0.5251, 0.5167, 0.4926, 0.4784 
0.4527, 0.4361, 0.4256, 0.4151 

(2,0! (0,2) 1-2,0 1 1-5.4) 

0.5469 

0.3969, 0.3711, 0.3595, 0.3424 

1Y1TV4.4) (-U) _(-4,4)_ 

0.4199 

0.3320, 0.3170, 0.2909, 0.2778 

1 i ,1) i i uN i I M f (~')A ) 

0.3061 

0.2390, 0.2118, 0.1801, 0.1643 

(0.0) 

0.0000 

0.0000 


4.4 Calculation of the two-magnon Raman intensity 

As before, [fij (also see Chapter 2) we can express the propagator in terms of it’s non- 
interacting components G 0 M- The mode assignments allows us to write GoM in a 2 x 2 
mat rix form with the mode energies and the symmetry factors labeled by the mode numbers 
,i( /,',TT',4 b Hence we can sum up the perturbation series in powers of the interaction 

between the magtions and which yields, 


[CM] = 


[Con] 

1 + J[Gon] 


(4.19) 


As in reference [(>] the above equation is written in terms of the matrix elements, 


[Gt)n(u0]-4.4 a n , 

[GonMJfis = Bn and 
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[GonM],4B = [Gon(w)]s<4 = C n 


0' n (u>) 


An±Bn + 2Cn + 2 J(A n £ n - Q) 

(1 + JA n )(l + JB n )~ PCI 


where ,4„. and C n are defined as, 


(4.20) 


.In 

/>’, 


1 ^ , , ( 2.1 
:;i^ rn °n ( — 

'< „ \U-’ n 


-f^Y 


2JJ 


a n u\ 1 / oo V 

IT ) + 2 V2j; 


i 


i 


1 v , , (2 J 

"■ - 2_ ™ °n ~ 

r > V W,i 


l V ,,72J^ 2 

~ 2- m 


ta 


n \ 2 7 flnW 


’i 2 2* 

[ b M 


2 JJ 
1 


+ 




2J J 2 \2 JJ 
1 


uj — 2u n t o + 2u n 

1 1 


La; — 2 tJ n u + 2 u n J 


La; 2 u) n uj + 2 lo u J 


where 


w -v 1 


> ‘J. 
-a* • 


{3 4- |cos(3 ; n I ) cos(^n y ) + 7 2 } and 


; , , , lid 

» 1 ~ Hr 


In Fit’;. (4.3) we have shown the Raman scattering spectra in the Bi g symmetry as a function 
of transferred photon frequency w (in units of J) for a system size 10 x 10 in presence of 
disorder (<7 ^ 0.1) and F// = 10. It is noticed that the linewidth to peak position ratio is even 
hh-her than unity, whereas the experimentally value is closed to 0.38[8]. This shows clearly 
th.i, the strength of disorder <r (or equivalently the mean square deviation {(t 1 ) 2 )) does not 

1 to be as large * 0.1 in order to get an agreement with the experimentally observed 

,]ma for the Human .pc hum. Fig. (4.4) shows the same plot for more moderate values of 
„ viz n - 0.01 and (1.0:1. The spectral lineshape for the pure case (<x = 0) is included for 
comparison. This plot is reminiscent of the magnon DOS shown in Fig. (3.4). There we have 
considered much larger strengths of disorder (a = 0.1 and 0.3). It can be observed that the 
linewidth to peak position ratio increases from 0.283 corresponding to the pure case to 0.37 
for ( T v 0.01 and 0.410 for a = 0.03 whereas the position of the peak roughly stays the same 
in nil die cases. So it looks like a - 0.01 roughly corresponds to the value that is needed 
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4.4 Calculation of the two 



I’ it’, lire -1.8: I he two-niagnon scattering intensity is plotted as a function of transferred photon 
frequency w' (in units of J) for the hopping disordered case with a gaussian disorder of width 

rr U.l. 

for ; ■ < : i ; i - ■ t r • : i . • m * \ disorder strength in order to obtain agreement with the experimentally 
obtained Hainan lineshape. Also it is seen that in presence of disorder there is a pronounced 
as\ muiet ry in the higher frequency region of the spectrum and a considerable intensity is 
observed all the wav uptu (i./ and even more. This can be explained by arguing that the high 
energy magnon inodes are strongly affected due to a cooperative effect arising from local 
cort e!a? ions of hopping disorder, which resulted in an appreciable one-magnon DOS beyond 
the mtidmtim energy ‘2J for the pure system[3]. So it seems that a much lesser value of 
the disorder strength can actually account for the puzzling features of the Raman lineshape. 
Finally, in Fig. ( 4 . 5 ) we plot, the scattering spectrum in the strong coupling limit (U/t = 100) 
for various disorder strengths viz. a - 0.01, 0.03 and 0.1 in order to draw comparisons with 
thi* results obtained bv Nori t‘t e/.(l] since their work deals with the strong coupling limit i.c. 
r/t ... x . The linewidth to peak position ratio in this case comes out to be 0.333, 0.406 
and 1 corresponding to rr = 0.01, 0.03 and 0.1. Also there is a pronounced asymmetry with 
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i ('spec t t o the two-magnon peak in the higher frequency regime. To obtain an idea about 
the asymmetry we measure the quantity /( u = 4J)// mar (u:). This has a value close to 0.4 
which agrees \eiy well with the corresponding quantity obtained from the experiments done 
on La-jCuOi (see 1* ig. (1.2)). Nori et al. obtained agreement with experiments for aj ~ 0.5J 
(they have defined crj in units of J) which is much more than what we have considered. 
This can he explained by noting that their Quantum Monte Carlo (QMC) calculations of the 
Hainan spectrum tor a N = 1/2 Heisenberg antiferromagnet in the B lg scattering geometry 
reveals a much lessor linewidth to peak-position ratio (~ 0.176) since in their calculation the 
iitlnif'ie mae.non damping has not been taken into account (remember we obtained 0.283 for 
l he corresponding quantity from Fig. 2.2 which is much closer to the experimental value 0.38 
than that of Nori 1 1 al.). This explains the reason why they need a much bigger width (cr) in 
the distribution of the exchange constant J so as to have an agreement with experiments. In 
fact, the evaluation of tin* different spectral moments of the Raman lineshape due to Singh 
,t shows the ratio of the second moment to that of the first one (which is like the 

hnewidth to peak-position ratio) as 0.23 which tallies very well with the experimental value 
II 3\ In Appendix B we have included a brief discussion on the evaluation of different spectral 
moment which was extensively used by Singh et a/.[8, 9] to analyze the Raman lineshape. 
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Fii-nn* 4.4: The two-magnon Raman lineshape is shown for more moderate values of disorder, 
\i/. a 0.01 and 0J1T Tin 1 pure case (a = 0.0) is included for comparison 


In conclusion, \v«* have studied the effects of hopping disorder on the two-magnon Raman 
interisit v using the rnagnon energies calculated numerically. The energies for the disordered 

assumed to have a one-to-one correspondence with the ones in the pure case. This 

correspondence is established by assigning mode numbers n x and n y which are related to 
,hr momentum labels via the relation q lh = (f)n s/y . The validity of the assigning scheme 
t becked for the pure case by re-evaluating the mode energies using the mode numbers 

;ui q -marine them with the ones obtained numerically. Once the consistency of the 

a-iiMiiue scheme is checked, the* numbers can be used to determine the symmetry factors 
which are required for calculating the Raman intensity in different scattering geometries. 
For U„. Hainan lineshape in the B„ in presence of a small amount of disorder (<r as small 
as |j ill i symmetry it is observed that the lineshape broadens and alongwith there appears 
a substantia! .. mttmity towards the higher frequency side while the low-energy modes are 
wenkk all'.-, ted. it is also observed that the spectral intensity persists even beyond (,J where 
tla- join, nMi'.noji DOS diverges at 4 J. We also obtained a good agreement for the linewidth 





®/J 

I-'ii’in *• t . : Tin* hV(Hiia,i’,nou Raman lineshape is shown for the strong coupling case (U/t = 
100i in older to compare with the results of Reference [1]. The values taken for disorder are 

a 0.111 , U.Od and 0. 1. 

and aammetry of the fwo-magnon profile observed experimentally from the Raman studies 
on cuprate insulators for <r - 0.01. T'his clearly reflects the role of zero-point fluctuations of 
the lattice to explain the anomalous features appearing in Raman spectrum for the cuprate 
ant iter mm, umefs via a highly asymmetric magnon energy renormalization arising due to 

correlation of hopping disorder. 
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4.4 Calculation of the two- m agnon Raman intensity 

Appendix A 

Herr we lind an expression for the symmetry factor appropriate to the Bi g scattering geometry. 
’I he Henry- London Hamiltonian for the B lg scattering geometry is written as (see Eqn. (4.5)), 

II = .4 S(r).{S(r + px) — S(r + py)} (4.21) 

r,p=±l 

I he two mutuum propagator in the non-interacting limit for this case is, 



whose r. r' denote lattice sites and r = r + 5, r' = r' + 5 ', where 8 / 8 ' connect to the sites 
neiidihourini'. to r and r' in all the four directions i.e. 5 / 5 ' can be px and py with p = ±1. 
Also the spin raisine and lowering operators forming the propagator is shown in the Fig. 4.5. 
Cleat!;. ‘ consists of Hi terms with 5 and 5 ' having all possible combinations which are 

enumerated below (the energy dependance is dropped temporarily). 


\H u (i\r') “ \ + '(r + j,r' + i) + x + '(r + i,r'-r) (4.23) 

- \ + "(r + £,r' + ij)-x + ~{r + x,r' -y) 

4- \ + "(r~ i,r' + i) +x +_ (r-^,r / -x) 

- \ + “(r ~i,r' + y) -x +_ (r-x,r'-y) 

- \ + ~{r + y,r' + x)-x + ~{r + y,r' ~x) 

+ \*~(r + y, r' + y) + x + "(r + y, r' - y) 

- ,\ +- (r-y,r' + x)-x + -(r -y,r'-x) 

+ x + "(r-y,r' + y) + x + '(r-y,r'-y) 

O obtain the symmetry factor one has to fourier transform all the 16 terms in G 0 (w). As an 

xeunple, the calculations for the first two terms are shown below, 
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Fin 


i i 't Tin* fin 


Tuktin’ the 


Which < .Ul 


( 2 1 Aik! tin 


Fourier 1 n. 






r + 5 


r’ + 5 




1 + + 

~ry v Sr + x “ Sr + y) 




o ( S r ’+ x S r’+ y ) 


nr S.ti: The spin raising and lowering operators are shown for the Go(H). 


t term is. 


+ (r,r'),\ + (r + x,r' + x) 

r,r' 


Fourier Transform, 

d 




r*r 


he written as. 


i J %Ex-+to)x*-l-i) 


second term. 


r/“ Ex^Cei-V-K + i.r'-i) 

^ j* t r / 


f. •• the above quantity yields, 

d. 


f y" y^ ^-+^) e i?.(r-r , )^+-(g') e »V-(r-r , -2x) 

I * J IW , mmmd Attmm u* 

J #T w, „* /»* 


r,r 


,*:'\:devt to. 


Which i:- 
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" 9 

Likewise we ran calculate tilt' contributions from the other 14 terms and all of them combine 

to t’i 


( \ + (9' w i)\ + -w 1 )[(cosg r - cosq y ) 2 ] (4.24) 

Ib'urc the sutusietry factor used in the Bi g geometry is (cos q x — cosq y ) 2 . By taking the 
appropriate sie.tis ( posit ivt*) for the spin operators corresponding to the lower line in Fig. 4.6 
for the A k. geometry. it is trivial to show that the symmetry factor is (cos q x + cos q y ) 2 . 
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Appendix B 

Smuli ft al.i s! '"dculafed the different spectral moments of the Raman scattering spectrum 
which are obtained in terms of the multiple spin correlations in the ground state of the 

I lei wnberg bang Hamiltonian, 


" = j, E Sf.-ff + ■!., BSfSf + SfSJ) (4.25) 

<n) (u) 

where the U limit corresponds to the Ising model and is trivially solvable. The authors 

have expanded the frequency moments of the Raman lineshape in powers of x(= pf-) and 
extrapolate the series to estimate the quantities of interest for the Heisenberg model (x = 
1 q't . In pai'tit uhir. the scattering intensity is given by, 

/!-') = E i<0|//R|i}| 2 <^ - {Ei - E 0 )) (4.26) 

v. he; e is ‘he '.'round state energy, E t is the energy for the state |i) and #r is the FL 
Haimhouian lie's ribed earlier. The spectral moments are extracted using, 

Mi = J ljI[u >) du (4.27) 

*il U L 

M/.r =/(-•-• W.r^p for n>l (4.28) 

•.viieje ! : i'. the urn-rated intensity given by l-x = f I(u)du. 

A comparisosj : . done for the B lK scattering spectrum in La 2 Cu0 4 with the various cal- 
culated m on:e,u!s. j ; w . [ )tM k frequency u> p for the two-magnon excitation process is somewhat 
smaller Hum the !«M moment .U,. The second moment M 2 measures the width of the spec- 
trum, hhe ruuo Mi comes out to be 0.23 which is good agreement with the experimental 
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Table 4.3: experimental values for the peak frequency, first, second and third moments of the 
/ 1’ 1; seatterim’, piotile in I,a 2 C'u0.i corresponding to the laser excitation wavelength of 4880 
,4 n . 1 he experimental lineshape is shown in Fig. (1.2). Also the values calculated using 

equations (4.271 and (4.28) are given. 


Mi 

Ah 

Ms 

M 2 /Mi 

Ms/Mj 

! cm 1 

107)0 cm -1 

1000 cm -1 

0.28 

0.27 

.6 ./ 

0.8 J 

1.0 J 

0.23 

0.27 


value ( -- 0.2 s ' u 1 he ratio \hf.\f i is the measure of skewness of the intensity profile which 
put hound/, on the hieh frequency tail of the spectrum. This ratio comes out to be 0.27 which 
, U! excellent experimental agreement. In the Table 4.3 we present the data from reference 
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Chapter 5 

Raman spectra in presence On-Site 
Disorder in Antiferromagnets 

5.1 Introduction 

In this chapter, we focus our attention on a class of correlated electron systems which are 
intrinsically disordered. The reason because of which they are of great importance is that 
they have tin* character of both Mott-transition and Anderson localization^]. In some cases, 
such as in U, x S:\YO.» and V a 0 ;{ , the metal-insulator transition is accompanied by a loss 

of ,\F order'd, :{j whereas in NiS 2 _*Se x one observes AF metallic phases [3]. 

Recently. interesting issues such as the interplay of quantum fluctuation and disorder ef- 
forts were studied] tj. In particular, questions like the following are addressed and answered, 


• Are transverse spin fluctuations quantum corrections to sublattice magnetization en- 
hanced by disorder? 

. Dues gapless AF state* exist? What is the critical disorder strength for which the AF 

lung range order ( AFLRO) is destroyed? 

More importantly for the* present discussion, the magnetic excitation is studied in presence 
disorder and found that the long wavelength, low energy magnons are weakly affected. So 
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those modes continue to be well defined excitations even in presence of disorder. However, 
it is observed that the high energy modes with energy ~ 2 J are damped. A self-consistent 
evaluation of magnon damping is necessary since the damping is found to be proportional to 
the maguon DOS. which diverges at the upper band edge 2 J. This leads to a self-consistent 
magnon damping which is proportional to the fourth power of the disorder strength and 
hence quite small. Since the two-magnon Raman scattering process probes short wavelength, 
high energy magnon modes, this magnon damping will be relevant in analyzing the Raman 
lineshape in presence of on-site disorder. 

The outline of the chapter is as follows. In the next section we mention the essential 
result s[-l] of tin* disorder induced corrections coming out of the perturbative analysis in the 
strong correlation limit. In section 5.3 we evaluate the self-consistent magnon damping which 
is obtained by considering the vertex correction which occurs at the second order level. This 
damping is later utilized to obtain the Raman lineshape for disordered antiferromagnets. 

5.2 On-site disorder 

In this section, we mention a few important results obtained by Singh et al.[ 4] via a perturb- 
ation calculation upto second order in the disorder strength. We write down the following 
Hubbard Hamiltonian with random on-site energies and with a filling one-fermion per site, 

H m £ c<n,> - 1 53 + b.c.) + U S n< t n 4 ( 5,:L ) 

«> (ij)o * 

The random on-site energies are chosen from a distribution where -W/2 < u < W/2, the 
distribution width W parameterizing the disorder strength. 

Treating V s £, an* in the Hamiltonian as perturbation they obtained the correction 
to the particle-hole propagator *, diagrammatical^ where the diagrams are shown in Fig. 
(5.1). In the first order they obtain, 
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fe] (U = 0 (5.2) 

owing to the particle- hole symmetry in an antiferromagnet. For the second order corrections 

the following results are found, 


r r i2 t 2 1 (e, £ i+cS) 

jj - 2 

[t i2 ^ ^ (5* £t+j) 

L^X'oJij _ ^34 JJ 2 


(5.3) 


Next the configuration averaged self-energy (J2) is obtained using (]£) — U 2 (6\ 0 ) where the 
diagonal and the off-diagonal terms can be written as, 


<S«> = -a (5-4) 


i„ units of ‘i'r and <7 being JgJ. This yields a momentum independent renormalization 
,,, the magnon propagator and hence to the energies (3 = J( 1 + c)) indicating an effective 

stiffening of the magnon modes due to disorder. 



5,2 On-site disorder 




c») (j) 00 w 

Fimiir "t . 1 : Th<- first (a ■ d) and second order (e - 1) diagrams shown for calculating zeroth- 

nrd<T partid** hole propagator [ \ o] - 
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5.2 On-site disorder 


T1 "‘ :i, " h " rs h “"- »!«• ‘■'iln.lated the first order correction to the eigenvalue A(u) due to 
,h ' 1 tV] "' hirh lhl ' knowledge of the eigensolutions {A™, |,)} of |\V)]. which has 

t r.-iti'-Kit tonal symmetry within the two-sublattice basis. They are given by. 


aiv> - 


U A 


jd - + 7, 2 ) 



— sin 0/2 
cos 0/2 





Vi- w 

x/TTU 


= *g-r, 


an<i I s th<* ’/<Toth ordrr particle-hole propagator given by, 


[ fr ’o(^)] 1 = 


1 -f- CJ 

7<? 


7g 

1 — u 


in enemy units where 2./ = — 1 and = ^(cos^ + cos<fo). 

So the first order correction can be obtained as, 


(5.5) 


(5.6) 


'WS ,, = « ) I(^ 0 ]I<>) = -<T(1-7, J ) 


(5.7) 


wheje c ‘ \ 1 .* . The magnon energy for the mode q is then given by the solution of 


l- v /^+ 7 2-aJ 1 > = 0 (5.8) 

l his readily yields. 


AjU = (1 + 2<r) v /l -% 2 = (1 + 2 a) u q (5.9) 

V. here we ha\e defined *•, ~. J\~l q - This tells us that the magnon energies get modified 
by a momentum independent multiplicative factor (l+2<r) and hence the modes get stiffened. 

Moie significantly, tin* authors have calculated the vertex corrections to the magnon 
y ■ , .* which results from a second order scattering process, yielding a magnon damping 

at mder ( ). 1 his was neglected in the earlier discussion where only configuration averaged 
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self eneig\ l\M) " a - s <■ onsidered. Our idea is to utilize this disorder-induced damping in the 
two magnon pi opngntoi and hence obtain the Raman scattering intensity. 

5.3 Magnon Damping 

Heir we will evaluate the damping of magnons resulting from disorder-induced scattering. 
For th, tt we consider the following configuration averaged self energy at the second order 

lesel.g; 


r = OJG'oS) - <S)Cr’ 0 (S) (5.10) 

who h i:n > n poraies vertex corrections. Next the magnon energy renormalization is considered 
tine to the self energy F upto order IF 4 . To this order, it is sufficient to examine the lowest 
oii'ei i orreelion (</jF|</) to the eigenvalue of U{ 1 — UXp ^ re ) matrix for the pure AF. 

So the eigenvalue correction comes out to be (see Eq. (14) of reference [11]), 




hAj/M) - « i) 


• cos i 


sO{l/i\)J2 M “MX 1 -7*)/(w* -w 2 ) 


2 / " 1 ( a. i — cr.i ) cos c/j. cos q y (l — cos k x cos k x ) / (u>l — u ? ) 
Z 1 1 (it.j - <r \ ) cos 2q x cos 2q y ( 1 - 72 A :)/ l ~ uj2 ) 


(5.11) 


whose / is the coordination number (=4 in 2D). The second term in the LHS which is 
suhti.u ted from the first one is the correction to the eigenvalue in the long wavelength ( q < 1) 
, it: d ;.,v. enerjw u .. i i (i.e. terms upto order q 2 and uj 2 are retained) <r 4 and cr 2 denote, 

o. X, and . T; W where 



(5.12) 


h * x tffMp >U* * f i ’< ’JAT j. I' 
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To obtain the magnon energy renormalization for a mode q we need to solve the following 

equal ion. 


1 - yo + g - «('l - «M(u,) = 0 (5.13) 

whore h M, < T l 1 "7 1 is the eigenvalue correction due to first order self energy (S) and it’s 

odor! is db-euv-ed earlier. Here, we only focus on the imaginary part of the magnon energy 
renormalization in order to obtain magnon damping. So we neglect the terms in the above 
equation whi< h contribute to the renormalization in the real part of the magnon energy. So 

the above equation reduces to. 


1 - \p + 3, ! - = o (5.14) 

when* lab vV.-'u-M' ■ • Solving this to first order in q q we obtain, 

.i, = «, - hi = u +.T, (5.15) 

' U> q 

v. here x ; 'l q i'U is taken to be 1) and r,(w) = is the magnon damping. This 

;e]; . us that the ruaenon damping P„(-;) can be obtained from the imaginary part of £A< 2) M- 
We V, ill diow later it! this section that the calculation of the imaginary part involves of r,(w) 
whit h is therefore ha> to he determined self-cons istently. So we look for a self-consistent 

uqiialinn fur I 
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Hero, we consider the imaginary part of <SA< 2 >( W ), 


I,n[(U 2 M] = ^<-os^(l/-V)I» 1 [L K -u; 2 ) 




(5.16) 


* ) 7 *■ l / 

“ ' \°d — ^2 ) cosq x cosq y 


(j_ — 7fc) 
K-uH) 


(cr-i — oj|) ^-(cos 2q x + cos 2q y ) 


(! - 72 *) 




This can he written as. 

ImTAp.o'p •- - a.' cosdl 1/A )I”i 


T. - ct 1) 


(1-7?) 


L k 


2 LJk l UJ - Uk 


i i— ) 

-Uk 0J + (Jju J 


U)+U)k) 


( 1 

1 i 


m + ui k J 

1 

i — i 

( 1 

2u/~ 

IcJ - LJje 


(5.1' 


+ Wjt 


Since .0 is positive the second term !/(«.' + Wfc ) does not contribute to the imaginary part. 

Hence we get , 




v.d(l/. V ) 




1 k 


2 oJk 


(5.18) 


I he 't ! it »!!•' appe.u iug in 

ha’. !!.»’ the ft .nil. 


2 Z" 1 ( a -i ~ orT) cos q x cos — 5(u> - u k ) 

ZuJk 

-Z~ l (<*4 - cr 2 ) i(cos2 q x + cos 2g y ) — ~ w *) 


the above equation' can be replaced by a lorentzian of width T g 


$(u/-w fc ) 


1 


(5.19) 


S nn t* ai ** ij/r: »*\t i*« ! nt 


TT (ul - Wfc) 2 + T 2 

, j fc- values which are close to the wavevectors <?, the choice of the 


width in Jim. case appears to 


a be appropriate. This allows us to write, 
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r = 

1 <i 




= cos«(i/iV)fx; K-o- 2 ) 


- 2 Z (a 4-0-3) cos q x 


2a;, (^-o;,) 2 + q 

(1-7D r, 


( 5 . 20 ) 


cos q y - 


2u k [u - Uk) 2 + n 


■ z “ + jVIWr p 


which is solved self-cousistently for T,. In Fig. (5.2) we have shown the contour plot of 
magnon damping vs q x and q y for W/U = 1 where W denotes the disorder strength. On the 
(</,., q v ) plane constant values of magnon damping are shown for one-fourth of the Magnetic 
Brillouin Zone (MUZ). It is seen that the damping is maximum along the MBZ edge precisely 
showing the constant value where the magnon energy is maximum (= 2 J). This tells us that 
the damping is maximum for the high energy modes which is expected since the damping is 
found to he proportional to the DOS which shows a divergence at the upper band edge (2 J). 
This damping is subsequently utilized to determine the scattering intensity in the strong 
coupling limit (I'/t = 100) for two different disorder strengths viz. W/U = 0.5 and 1 (Fig. 
5.31. It is observed that the lineshape broadens. In view of the divergent magnon DOS the 
damping obtained here is substantial even though it is proportional to (cr 4 — cr 2 ) which comes 
from vertex corrections and this results in an observable disorder-induced broadening to the 
scattering intensity. It can also be observed from Fig. 5.2 that the low energy modes are 
negligibly renormalized and hence there is an insignificant change in the low energy regime 
of the two- magnon Raman spectrum (Fig. 5 . 3 ). 






5 3 Magnon Damping 


108 


In conclusion, wo have studied the effect of disorder-induced magnon damping on the 
two-magnon Hainan lineshape in an antiferromagnet. For the high energy modes (close to 
2,1 \ which are particularly relevant for Raman studies, the damping is proportional to the 
piodm t of DOS and a term which goes as fourth power of the disorder strength. Because 
of the divergence of the magnon spectrum at the upper band edge these modes are affected 
the most. This shows up in the two-magnon scattering intensity where there is a noticeable 
renormaliration in the spectral intensity towards the higher frequency regime. 
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Chapter 6 


The Real-Space Formalism 

/ 

0.1 Introduction 

h.uiici we haw calculated the Raman scattering intensity corresponding to the pure case in 
( h. spier ... i here we described in detail the perturbation theoretic scheme for systematically 
iTt, I!!!!!!!’, the two magnon Raman propagator and hence the scattering intensity from it. It 
w.f, id o shown that the t wo-magnon propagator in the non-interacting limit can be expressed 
in terms of two transverse spin propagators [x~ + ] and [x + ~] for which the translational in- 
wiriunee holds within each sublattice and can be represented in terms of two 2x2 matrices 
; \ * I, </. ~')j ami ( \ + ~ ( — u>i )][1 , 2]. The calculation of the Raman intensity requires the 

knowledge of the spin- wave energies which are obtained using the relation u q = DJ^Jl — 7^ 
where 1) denotes dimensionality and 7 , = J2tLi cos q v /D. In presence of hopping disorder 
M 'h.tpfer •{) again we have obtained the scattering intensity by using a scheme of mode as- 
"iimments. This essentially establishes a one-to-one correspondence between the pure and 
t !a* disordered system. In order to establish the correspondence we calculate the spin-wave 
energies using the exact-eigenstate analysis in presence of disorder and assume that these 
exact, energies can be labeled bv the mode numbers n x and n y which are related to the mo- 
mentum labels q r and q y via the relation n x / y = ~ q x / y where L denotes the size of the system. 
Now this kind of an attempt to establish a correspondence will be appropriate only if the 
disorder strength is weak since we will still be assigning the same momentum labels as the 
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pun' system. This assignment scheme will not work if the disorder strength is moderately 
high. 


In this Chapter we will formulate the problem of obtaining the Raman intensity entirely in 
real space i.c. using the eigensolutions of the zeroth order particle-hole propagator determined 
from the numerical approach which is discussed earlier in chapter 3. The construction of the 
zeroth-order two-magnon propagator [G°(tu)] in terms of eigensolutions of [x°] is described 
in section 6.2. Here we adopt the method described recently[3] for obtaining transverse spin 
correlations in terms of magnons amplitude and energy. The method is illustrated in Appendix 
A following reference [3]. In section 3 we discuss the numerical method of calculating the 
Raman scattering intensity which is obtained from the imaginary part of the two-magnon 
propagator. Also we mention that the calculation of the scattering intensity is numerically 
intensive and this grows considerably with system size. 


C.2 The two-magnon propagator in Real Space 

It is shown earlier in chapter 2 and reference [4] that the two-magnon propagator can be 
obtained from the correlation function (pj Pj,) where P- 5 = E r S( r ).S(r + 6). In the non- 
interacting limit., the propagator can be written in the real space as, 


[C"M] = i / ^-(x- + (r, r>i)][x + -(r + *.*' + s ’’ a ~ “">! t 6 ' 1 ) 

Where t.l«. indices (r.i) (or (i'J')) denotes that there are 4 choices of S (or S') via. ±x,±y 
f„r a particular site r (or r'). For convenience we define new indices „ and yl given by 
/i & (r.ij ami // = (Cj') and from now on we will denote the the propagator as [G°(w)]„,/. 
Corresponding to a particular location denoted by the indices (r.r'J in the [G MW matrix 
is a 4 x 4 matrix for all possible combination of S and S'. This tells us that the dimension 
of the matrix [GV)1 « g^n by [r® S[ x [r'O S’] where r/C runs over the number of lattice 
sites ami US' connects to all the nearest neighbours and hence we have a 4iV x 4iV matrix. 

Kq. (6.1 ) can also be written as, 
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6.2 The two-m 


[G MW - ij ^[x + ( r ? I * , ) a; i)][x + (r + J,r' + 6',iOi -to)] (6.2) 

\\ Ik u [\ ] s aie the magnon propagators at the RPA level and are represented as, 


[\" + M] = 


[xoM] 
i - U[ x oM] 



-Mcu) 

1 - U\ n (u>) 




( 6 . 3 ) 


\\ hem { A„ , 6 n } are the eigensolutions of [xo(cv)] which are obtained from the exact eigenstate 
analysis and n denotes the labels of the mode energies obtained by solving the equation 
1 — l A„(u.’) = 0. Putting these back in the expression for the non-interacting two-magnon 

propagator one obtains, 


[( 7 °( uj )] 



]<Mr))(<Mr01 Y . |^ m (r + ^))(^(r' + ^)| 
Uil-UXnfa))*? Uil-UXnfa)) 


( 6 . 4 ) 


Where Ay = ay — u. Now the frequency integral has to be done to evaluate the propagator. 
Kach of the terms above will have two sets of poles corresponding to the advanced and re- 
tarded part of [\~ + ]. For convergence the retarded part has poles below and the advanced 
part, has poles above the real-ce axis. In the first term the poles are at ay = +u n — irj and 
d it) and for the second term they are at chi = +u; m — it) and — uj m + irj. So to have a 
non mro contribution in the frequency integral of the above equation it is necessary to include 
one pole above and the other below the real axis. Thus there will be two terms coming out of 
the integral, one consisting of the combination u> n - irj and -u> m + irj and the other — u n + irj 
and Uv, - it)- The first term can be seen explicitly if we Taylor expand the eigenvalues 
A,, (ay) and A m (uy) around the pole frequencies uy = u n and uy = -u> m respectively. This 
gives. 


. f d tl ^ l<ftn(r))(0w(r')| y- l<Mr + ftKAnfr + *)1 

1 J 2rr V -U 2 { uy - ce„ + irj) (^) wi=Un m -U*{& 1 + ~ irj) (§r)* 1= -, 


( dearly the above terms have poles at uy = u; n - irj and uy - a? m + irj. The quantities 



6.2 The two-magnon propagator in Real Space 


113 


and (^f) 


w 1= - Wm are s ^°P es °f £/(l — U A„(o; 1 )) vs. u;i evaluated at u; n and 
f ( 1 f A, n (a. 1 ! )) vs. d’i evaluated at — u> m respectively. The calculation of these slopes are dis- 
1 u,SSt ^ ( ^ a ^ 01 ately in Appendix A. To make the notations compact we define —U 2 

i,s ’ a " d - (;2 (^L=-„ as AW. 


U/ri 


Now the frequency integral is done by including the pole at lo\ = u) n — irj which gives, 

_ l<Mr))(<?nL(r )| y [0m(r + + 5')\u n -w 

n An * m Am^(tU — Ul n — -f * 77 ) 

where is approximated as irj. In order to get the second term in the integral again we 
Iaylor expand the eigenvalues A„(cj 1 ) and A m (a> 1 ) around the poles ui\ = —u n and clq = u m 

which gives, 

y- l^n(r))(0 n (rQ| y l0 m (r + l))((?!> m (r / + (5 / )| 

As before we can define — U 2 _ as A^ and — U 2 . _ as A^h Finally the 

frequency integral is performed by including the pole at uq = — u> n + irj which yields, 

y lPn(r)){^nl- u >„(r / )l y l^m(r + QKM r ' + 

n Al 2) m Am (W — U) n — U m + if)) 

rombining the two contributions, the two-magnon propagator at the non-interacting level is 

written as, 


rWW , , [ yp l^n(r))(</>n(r / )l yv 1 4> m(j + ^))(0m(r + £ )1 5' 

= i ? ZF A aS>(. - „„ - ^ + in) 1 

y IjnWHgl y ]^m(r + ^))(^m(r / + ^ / )l ‘ 

n An^ in A$(w + cu n + u m - irj) - 

Now the full two-magnon propagator is calculated by summing up a perturbation series in- 
volving repeated magnon interactions which yields, 



6.3 Numerical approach 


115 


6.3 Numerical approach 

The formalism described here is useful in calculating the two-magnon Raman scattering in- 
tensity for systems in presence of disorder such as, hopping disorder, potential disorder, 
impurities, vacancies etc. Since the method involves calculation of the scattering intensity 
entirely in terms of the real space representation of the propagators, it is applicable to any 
strength of disorder and a whole range of impurity (or vacancy) concentration. 

Here, we discuss the different steps of the numerical exercise that has been followed in 

the calculation of Raman intensity. 

• As a first step the Hartree-Fock (HF) Hamiltonian can be diagonalized in presence of 
disorder to obtain the eigensolutions {A n ,<£„}. 


• The eigensolutions are to be used to construct the (xo( w i)] matrix which should then 
be diagonalized. The eigenvalues of [xo] are linearly interpolated to calculate the mode 
energies, u>„. At the same time the slopes A^ and A^ can be stored. 


0 In a similar manner the matrix [xo(^i)] can t> e formed and the mode energies u m are 
calculated. Also we can store the slopes A^ and A^h 


. Having found the mode energies, [xo] matrix is to be reconstructed at the frequencies 
wV ^ and -w n - w which is subsequently diagonalized to evaluate the ei- 

genvectors |0„(r)) and |^(r + J)) at these frequencies. 

. If, '’MU is formed according to the equation (6.5). It is to be noted that [G»W 
i, » matrix of dimension 4 N x 4JV since we have four choices each for f and S' cones- 


6.3 Numerical approach 


116 


ponding to each r and r'. 


• I he matrix [G°(u;)] PiM i is inverted calculated using a standard NAG (Numerical Al- 
gorithm Group) library routine. 


• In the next step the inverse of [[G°(u>)] 1 + J 1] is to be calculated in order to obtain 
the full t,\vo-magnon propagator G(w) (= [[cvr +n]-' >. 


• In order to calculate the scattering intensity we need to do the following. For a particular 

choice of 8 and 8' consider the sum over all the site indices r and r' in This 

sum is to be evaluated for all choices of 8 and S'. This will yield 16 terms. Taking 
appropriate combination of these terms and finally considering the imaginary part yields 
the intensity for different scattering geometries. 

• Finally one needs to calculate the average intensity taken over several configurations. 

From the above discussion it is clear that the whole calculation requires a considerable 
amount of computing time, since the construction of the [xo] matrices for several frequencies is 
repeated quite a number of times. Also the inversion of the 4 N x 4 N complex matrices which 
are done using standard library routine takes substantially long time. To add to that, the 
whole procedure has to be repeated for several frequencies (transferred photon frequencies). 
Finally, taking configuration average would imply repeating the whole procedure all over 

again for several configuration. 
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Appendix A 


I he spin correlations are defined as, 


(Sr(t)SJ-(O) = -;/ 




(6.8) 


Writing the R11S in terms of the eigensolutions of [xo(^)j, 


(s-( () s J v))=-.7|i: I 


A«M 


;4>n (u)<P n {u)e 


- iu>{t-t ') 


(6.9) 


U A n (w) 

Now for calculating the equal-time correlation function (t = i') we Taylor expand A n (u?) 
around the mode energies u„ and substitute in the above expression which yields for the 

MIS. 


■ t du y, .-.W(t-t') (6J0) 

( ! > j i s defined as the slope of A n vs. u> graph evaluated at uj = ce n . Now the integral is 
obtained by evaluating the residues at the poles. This yields, IX! 


, ^n^n^Q^n) (6.11) 

Where the eigenvectors |<£„(u> n )) are evaluated by re-constructing the [xoM] matrices at the 
mode energies ^ fl anti subsequently diagonalizing them. We now describe the procedure of 
obtaining the slop.' from our numerical exercise where one has to evaluate the sum, 


E 


d n (u ~ U> n ) 


( 6 . 12 ) 


Also (/. 


refers to the denominator 17(1 - UKM)- So we write the above expression as, 
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E 


I his ran he written as, 


£ 


d n (oj^u n ) 

(oj—ui n ) 


dn (u> n +duj) — dn (u/ n ) 

duj 


(6.13) 


(6.14) 


N>>w (/».(«;„) = 0 .since this condition 


gives the mode energies. So the above expression reduces 


to , 


E 


WK))(I<V») I 

dn {un-hdu/) 
dui 


(6.15) 


Defining the quantity in the denominator as A n which represents the slope of d n (u) vs. w at 


l^nK))(^K)l ^ 61g ^ 

n 

We ran calculate the slope while obtaining the mode energies u n . Since uj u is calculated from 
the n*ut of the equation, l — UX n (u} n ) = 0 which is solved by obtaining A n (o;) for several close 
values of a.* on both sides of the root and then linearly interpolating between them. Thus if 
i A i and (A*,uj*) arc two sets of values for two energies very close to and on either side 

of the root then u is determined from (see Chapter 3), 


A(w„) 


]_ 

U 


'V + — — rr( w " 


• UJ, 




■ (Jj, 


(1 ) 

nJ 


n 


(6.17) 


Here we can calculate the quantity which is defined as the slope A n (see Eq. 6.16). 

Comparing Kq. (6.16) and (6.11) we find that A n = ±U 2 ( i t) u=ojn which we can evaluate. 
This essentially does the frequency integral appearing in equation (6.8). 
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